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NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


MEMORANDUM 1-16-59A 


A NUMERICAL METHOD FOR CALCULATING THE WAVE DRAG OF 
A CONFIGURATION FROM THE SECOND DERIVATIVE OF 
THE AREA DISTRIBUTION OF A SERIES OF 
EQUIVALENT BODIES OF REVOLUTION* 

By Lionel L. Levy, Jr., and Kenneth K. Yoshikawa 


SUMMARY 


A method based on linearized and slender-body theories, which is 
easily adapted to electronic-machine computing equipment, is developed 
for calculating the zero-lift wave drag of single- and multiple -component 
configurations from a knowledge of the second derivative of the area dis- 
tribution of a series of equivalent bodies of revolution. The accuracy 
and computational time required of the method to calculate zero-lift wave 
drag is evaluated relative to another numerical method which employs the 
Tchebichef form of harmonic analysis of the area distribution of a series 
of equivalent bodies of revolution. 

The results of the evaluation indicate that the total zero -lift wave 
drag of a multiple -component configuration can generally be calculated 
most accurately as the sum of the zero -lift wave drag of each component 
alone plus the zero -lift interference wave drag between all pairs of 
components. The accuracy and computational time required of both methods 
to calculate total zero -lift wave drag at supersonic Mach numbers is 
comparable for airplane -type configurations. For systems of bodies of 
revolution both methods yield similar results with comparable accuracy; 
however, the present method only requires up to 60 percent of the comput- 
ing time required of the harmonic -analysis method for two bodies of 
revolution and less time for a larger number of bodies. 


INTRODUCTION 


As shown in reference 1, the linearized theory value of the wave drag 
of any arbitrary configuration in a steady supersonic flow can be calcu- 
lated from a knowledge of the second derivative of the area distribution 
of a series of equivalent bodies of revolution obtained from areas and 
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forces in oblique planes. For slender nonlifting configurations the 
wave-drag contribution of the force term is zero or is sufficiently small 
that it can be neglected. The zero-lift wave drag of such configurations 
is given to a good approximation by the supersonic area rule developed in 
reference 2. The infinite series form of the supersonic area rule in ref- 
erence 2 utilizes a knowledge of the first derivative of the area distri- 
bution of the series of equivalent bodies of revolution obtained from only 
the areas in oblique planes. In reference 3 the supersonic area rule was 
expressed in a form which permits a numerical solution of the infinite 
series from a knowledge of the area distribution of the series of 
equivalent bodies of revolution. 

In this report a method applicable to slender nonlifting configura- 
tions is developed which utilizes a knowledge of the second derivative of 
the area distribution of the equivalent bodies of revolution. The method 
is developed from the integral form of the wave-drag equation given in 
reference 1 for nonlifting configurations. By a simple approximation of 
the second derivative of the area distribution of the equivalent bodies 
of revolution, the zero-lift wave drag is expressed in a finite series 
which can be evaluated numerically. The finite series form of the present 
method is particularly attractive in its application to the calculation 
of zero -lift interference wave drag between any pair of components of a 
configuration. 

The usefulness of the present method depends, of course, upon the 
availability of methods for finding the second derivative of the area 
distribution of equivalent bodies of revolution. Reference 4 contains 
information from which it is possible to calculate the second derivative 
of the area distribution of bodies of revolution. A method for finding 
the second derivative of the area distribution (as well as the first 
derivative and the area distribution) of wing and tail-surface components 
of the equivalent bodies of revolution is presented in an appendix. 

The present method of zero-lift wave -drag calculation will be 
evaluated relative to the method of reference 3 with regard to accuracy 
and computational time required. Numerical values of zero-lift wave-drag 
solutions for several simple analytical shapes computed by the present 
method will be compared with numerical solutions obtained by the method 
of reference 3 an & with analytical solutions. Equations of the area dis- 
tribution, the first derivative of the area distribution, and the second 
derivative of the area distribution of two families ::f viaos or - 
surfaces will also be presented in an appendix. 
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SYMBOLS 


a ij 


A n (p,e) 

b 

bo 

c o 

c(t]) 


D(P) 

D(P,0) 

I 

l( a i j } 3 d -) 
K 

l( 3,9) 

l x (P,9) 

L(x,p.0) 


distance between the ith and jth values of coordinates in 
the free-stream direction (see eq. (4l)) 

coefficients of a Fourier sine series expansion of S' (x, [3,0) 

span of a wing or tail surface 

span of an extended wing or tail surface (see sketch (c)) 

wing or tail-surface chord at the vertical plane of symmetry 

spanwise variation of the dimensionless half-chord of a wing 
or tail surface 

lateral distance between the longitudinal axes of a pair of 
bodies of revolution 

wave drag or zero -lift wave drag of a configuration 

wave drag or zero -lift wave drag of an equivalent body of 
revolution of a configuration 

total number of points specified on a given curve 

integral function used to evaluate the zero-lift interference 
wave drag between a pair of bodies of revolution (see 
eqs. (40) and (4-3)) 

tangent of the sweep angle of the 50 -percent chord line of 
any "sheared" panel of an extended wing or tail surface in 
dimensionless coordinates 

length of an equivalent body of revolution in the positive 
x direction, l(p,0)+l 1 (p,0) = total equivalent -body length 

length of an equivalent body of revolution in the negative 
x direction 

l i ontal projection nic- .L LJTt distribution intercepted on a 
riven configure ■ y set of parallel ■ ■ -ique planes 

ni-'ent to ‘l"' 


y i V;‘r “ ■ T> Hi 11 i “l til 

i'rec -i.tre an; ciynH 
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s 

s(x,p,e) 

as"( X i ) 

s m"U) 

t O 

t(x,y) 

a(|,ri,K) 

3 

6 

A 

r(a) 

t(S,tj,K) 

9 

* 


reference area upon which drag coefficient or dimensionless 
area is based 

frontal projection of the area distribution intercepted on a 
given configuration by a set of parallel oblique planes 
tangent to the Mach cones 

increment between values of the second derivative of the area 
distribution of a configuration at points x-^.-l and 

mean value of the second derivative of the area distribution 
of a configuration (see eq # (45)) 

maximum thickness of wing or tail surface at the vertical 
plane of symmetry 

thickness distribution of a wing or tail surface 

Cartesian coordinates in the free -stream, spanwise, and 
thickness directions, respectively 

transformation defined by equation (C 17 ) which defines any 
point in the dimensionless extended plan form of a "sheared 1 * 
wing or tail surface as the local dimensionless chord 
station measured from the 50 -percent chord 

s/M 2 -1 

longitudinal distance between the lateral axes of a pair of 
bodies of revolution 

taper ratio of a wing or tail surface; ratio of the tip chord 
to that at the vertical plane of symmetry 

dimensionless Cartesian coordinates in the free -stream, 
spanwise, and thickness directions, respectively 

dimensionless thickness distribution of the local airfoil 
section 

dimensionless thickness distribution of a "sheared 1 ' wing or 
tail surface 

angle defining the orientation of the parallel oblique planes 
tangent to the Mach cones (see sketch (a)) 

spanwise variation of the dimensionless thickness along lines 
of constant percent chord 

tan’ 1 ^ cos G 
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Subscripts 


B body 

E, H body E or body H 

H(E) quantities measured in the coordinate system of body H but 

related by a transformation to the coordinate system of 
body E 

i , j ith and jth points, respectively, of the total number of 

points specified on a curve 

Z lover limit of integration 

n nth term of a Fourier sine series expansion of S ! (x,3,0) 

R reference 

t tip of a wing or tail surface 

T total 

u upper limit of integration 

W ving 

Superscripts 

1 differentiation with respect to a coordinate in the free- 

stream direction 

dimensionless value of any symbol not specifically listed as 
such (also see appendix A) 


METHOD 

Review of Basic Theory 


The lineari sod-theory value cf the wave drag of any arbitrary 
configuration in a steady supersonic air flow is stated in reference 1 
to be the average of the wave drag of a series of equivalent bodies of 
revolution and cm; be expressed analytically by 



where 


“(e- 0 ) ’ - hj 


i(P,e) oZ(P,e) 


-i x { P,e) a -i x (p,0) 


s"( Xl ,p,0) - A. L'( Xl ,p,0) 


S"(x 2 ,3,0) - 


JL 

2q 


L*(x 2 ,p,0) 


Znlx-L-x^ 


dx 1 dx 2 


( 2 ) 



The quantity S(x,(3,0) is the frontal 
projection of the streamwise area 
distribution intercepted on the given 
configuration by parallel oblique 
planes tangent to the Mach cones . 
Similarly, L(x,f3,0) is the frontal 
projection of the lift distribution 
over the given configuration measured 
in the same oblique planes . The 
primes denote differentiation with 
respect to x. The coordinate system, 
angles, Mach planes, and the direction 
of the net force are defined in 
sketch (a). 


In this report only slender nonlifting configurations will be 
considered. For these configurations L(x,3,0) can be neglected and the 
zero-lift wave drag of the separate equivalent bodies of revolution for 
such configurations is given by 

pl{ 3,e) pi(p,0) 

D(P>0) = - ~ / / S"(x 1 ,p,0)S"(x 2 ,p,0)Zn|x 1 -x 2 |dx 1 dx 2 

2jtJ -Zx(3,0) 

(3) 

The analysis in reference 1 imposes no restriction on the slope of the 
area distribution at the ends of the equivalent bodies of revolution or 
of the given configuration. It should be noted that this freedom from 
end-condition restrictions requires rigorous adherence to the concept of 
the oblique planes when determining the area distribution of the equivalent 
bodies. However, for nonlifting slender pointed bodies of revolution, the 
normal cross-sectional area distribution can be used if the end conditions 
of the bodies are such that 
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S’(-lx,3,0) = S'(z,3,0) = o (4) 

For isolated nonlifting slender pointed bodies of revolution for which 
the slope of the area distribution at the base is not zero, the normal 
cross-sectional area distribution can be used only if the zero-lift wave 
drag is calculated by slender-body theory. The zero-lift wave-drag equa- 
tion given by slender-body theory for these isolated bodies of revolution 
is given in reference 5 and can written 



where R(z) is the radius at the base of the body. The lower limits of 
integration, — Zi, imply merely that the origin of the coordinate system 
need not be chosen at the half -body length. For bodies of revolution with 
zero slope of the area distribution at the base, S*(z) = 0, the second 
term of equation ( 5 ) vanishes and the resulting equation is identical to 
the equation (3) for equivalent bodies of revolution, once values of (3 
and Q are specified. 

The area distribution of a single body of revolution is independent 
of Q (but may change with ( 3 ) and the linearized-theory value of the 
zero-lift wave drag is given by equation ( 3 ) without recourse to equa- 
tion (l). Furthermore, if the body is sufficiently slender that neglect- 
ing the change in the area distribution with (3 produces only a negligible 
effect on the zero-lift wave drag, the linearized-theory value of the 
zero-lift wave drag of a slender body of revolution with S ! (z) =0 is 
identical to that given by slender-body theory. The body fineness ratio 
required in order to justify neglecting the effect of 3 on the body 
area distribution has been discussed in reference 4. 


Approximations 


Aside from the approximations basic to linearized and slender-body 
theories, only one approximation is necessary in order to express equa- 
tions ( 3 ) and (5) in a finite series. This is an approximation of the 
second derivative of area distribution, S'^x ), 1 of the equivalent bodies 

of revolution. _ 

ifhe variables 3 and 0 have been intentionally omitted from the 
notation in the interest of simplicity and will be omitted throughout this 
report except in cases where these variables are necessary to the analysis. 
For bodies of revolution S"(x) = S"(x,p). If the body is closed at the 
rear, l is constant, and if the body has a finite base, l = l(p). For 
wings and tail surfaces S"(x) = S"( x >ft>6) and Z = z(p,Q). 
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Let S TT (x) be approximated by constant values over small equal 
intervals; that is, a smooth curve is represented by a broken curve con- 
sisting of adjacent equal -length steps as shown in figure 1. The value 
of the second derivative of the area distribution at a given station, x, 
is considered to be approximately equal to the value at the beginning of 
the interval, station x^, in which station x lies. This approximation 
can be expressed analytically as 

S tT (x) « S tr (xi) for x-^ < x < x±+ x 

If the difference between the values of the second derivative of the 
area distribution at two adjacent intervals is defined at the beginning 
of the right-hand interval as 

A3"( Xi ) = S”(x i )-S"(x i _ 1 ) 


then 


S M (x) can be approximated by 
i(x) 

S"(x) * £ £S"(xi) 


1 = 1 


for x± < x < x±+ x 


( 6 ) 


where the upper limit of summation defines the beginning of the interval 
in which station x lies. It should be noted that I points are employed 
to define 1-1 intervals for the total length (see fig. l). 


Development of the Method 


The present method for calculating zero -lift wave drag is developed 
from analyses which evaluate equations (3) and ( 5 ) in light of the approxi- 
mation of S"(x) described above. 2 All analyses in this section of the 
report will be made in physical or dimensional coordinates. Experience 
has shown, however, that in performing the actual drag calculations it is 
more convenient to use dimensionless coordinates. For this reason, pro- 
cedures are given in appendix A for making various types of configurations 
dimensionless . 

Single bodies of revolution, S ! (-Zi) = 0, S r (l) ^ 0.- The general 
expression for the zero-lift wave drag of a slender pointed body of revo- 
lution is given by equation ( 5 ). It is convenient to write equation (5) 
as 

^During the preparation of the present report a numerical method for 
calculating the zero-lift wave drag of bodies of revolution by electronic 
machines or desk calculators was published in reference 6. Application of 
the method of this reference also requires a knowledge of S"(x) (exact or 
approximate) and may also be used to evaluate equations (3) and (5). 



9 


D = Di+D^ 


where 


l pZ 


- • 27. 


S" (xx )S" ( x 2 ) Zn | X!-X 2 1 dXi&Xa 


( 7 ) 


■h-h 


and 


D 2 = q 


[S’(Z)]‘ 


Zn 


ihil J S"(x 2 )Zn|z-x 2 |dx 2 i 

-z 1 ' 


v 2n 3R(Z) 

First consider the expression for D]_. Equation (7) can be written 

Z 

•'-Zi 

where 

Z 


( 8 ) 


Di = - ~ j' S"(x 2 )F(x 2 )dx 2 


(9) 


F(x 2 ) = J S"(xi)Zn|x 1 -x 2 |dx 1 
-Z 1 


( 10 ) 


Using the form of equation (6) for the variable x x and using the fact 
that S"(x) is approximated by a constant over each interval of the 
distribution (see fig. l), one can write equation (10) as 


F( 


and equation ( 9 ) becomes 

I 


x 2 ) « y ^( xi ) r Zn jx x -x 2 |dx 

L — 1 , — v-t 


1=1 


x i= x i 


Z p Z p z 

A3"(xi) / S" (x 2 ) / Zn|xi-x 2 |dxidx 2 (ll) 

i =1 -2 1 x i =x i 

Similarly, using the form of equation (6) for the variable x 2 , one can 
write equation (ll) as 


I I 


Di * - y yAS’^Xj^^'^Xj) J J Zn |xx~x 2 | dx x dx 2 (12 


i=i j=i 


x 2= x j x l= x i 
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After the double integration, equation (12) becomes 


I I 


Di ~ - 


JL 

4n 


^ ^^"(x i )Z^"(x j )l 0 (x i ,x J ) 


( 13 ) 


i=i j=i 


where 


I o( x i^ x j) = U-Xq) 2 ^2n|z-Xi| - + (Z-xj)*" (zn|z-xj 

( x i" x j) 2 ( in l x i- x jl - |) 


(14) 


The order of the summation in equation (13) can "be interchanged. This 
permits parts of equation (14) to he evaluated as constants. Hence, if 
the contribution to equation ( 13 ) of the first term in equation (l4) is 
summed over i, the contribution of the second term of equation (l4) is 
summed over j and the contribution of the second part of the third term 
of equation (l4) is summed over both i and j, equation ( 13 ) can be written 


4n 


D X * - 


where 


and 


I I 


K i ^ AS"(xj )+K 2 ^Z^"(x i )+K 3 - ^ ^i5sS"(xi)/!iS"(xj )l 1 (x 1 ,xj ) 

(15) 


j=i 


1 = 1 


1=1 J=1 


Kn 


1 

= ^AS"(x i )(Z-x i ) 2 (^Znlz-Xil - ^ 

i=i 


1 

K 2 = ^AS"(xj)(Z-Xj) 2 (^Zn|z-Xj I - f) 

j=i 


I I 


Ka - | ^ ^AS"(x i )Z^"(x J )(x i -x J )' 


1=1 J=i 




J 


> ( 16 ) 


!i( x i, x j) = (x 1 -Xj) 2 Zn|x i -xj 


(17) 


Now consider the expression for D 2 given by equation (8). Using 
the form of equation (6) for the variable x 2 and the fact that S"(x) 
is approximated by a constant over each interval of the distribution, 
one can write the second term of equation (8) as 
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l 1 l 

G(x 2 ) = PP f S"(x 2 )Zn|z-x 2 |dx 2 * ^lli Y AS"(xj) f Zn|z-x 2 |dx 2 

^ 3 . J-l X J Me 


-u 

After the integration, equation (18) "becomes 

I 


g(x 2 ) * Sliil I AS" ( xj )J" ( Z -xj ) ( Zn | Z -xj | -1 ) 


( 18 ) 


(19) 


l=i 


Again, since S"(x) is approximated "by a constant over each interval of 
the distribution, it can be shown that 


-L 

^A£"(xj)(Z-Xj) » S' ( Z ) 


( 20 ) 


1=1 


and equation ( 19 ) becomes 


G(x 2 ) » pp ^AS"(xj)j(xj) - hhlll 


J=i 


where 


J(xj) = ( Z -xj )Zn| Z -Xj 


( 21 ) 


( 22 ) 


Substituting equation (21) in equation (8) 


D 2 


9 Jl 


it 2 


i [S'(Z)] ; 


Zn 


PR(Z) 


+ S 


1 

«(z)^AS"(xj)j(xj)j (23) 


J=i 


Finally, the zero-lift wave drag of a slender pointed body of revolution 
with finite base area, such that S'(z) / 0, is given by the sum of 
equations ( 15 ) and ( 23 ) as 


D* - f 

lit 


I I 


K x ^AS"(xj ) + K 2 ^AS"(x i)+ K3-X ^ AS"(x i )AS"(xj )l 1 (x 1 ,xj 


L- 1=1 


1=1 


1=1 J=1 


q 

Tt 



[S*(Z )] 2 



2 

m(i) 


i 



( 21 ) 


where K x , K 2 , K 3 , and I x (x x ,xj) are given by equations (l6) and (17)> 
respectively, and J(xj) is given by equation (22). 



12 


Single bodies of revolution, S'f-Zi) = S y (z) = 0 .- For a slender 
body of revolution pointed at both ends or with finite base area such that 
S y (z) =0 the last two terms of equation (24) become identically zero and 
the zero -lift wave drag of such bodies of revolution is given by equa- 
tion (15). Equation (15) can be simplified for S 1 ( Z ) =0. For this 
condition, it can be seen from figure 1 (since S"(x) is constant over 
each interval, S T (x) is linear over each interval) that 


I 

i=i 


I 

£>■(*,) = ° 

j=i 


for S T (z) = 0 


( 25 ) 


Hence, the first two terms of equation (15) vanish. Further, upon 
substitution of equation ( 25 ) in equation (20) it is easily determined 


that 


I I 

2^^6"(x i )x i = ) XJ =0 for S’(Z) =0 

i=i j=i 


( 26 ) 


With equations (25) and (26) it can be further determined from equa- 
tion (l6) that K 3 is identically zero. Hence, for S ! (z) = 0, 
equation ( 15 ) becomes 


Di ^ ^ Z Z^ ,,(xi) ^ n(x J )ii(xi ' x J ) (2?) 

i=i j=i 

Finally it should be noted that I 1 (x^,xj) = 0 for i = j (see eq. (17)) 
and in view of symmetry of the matrix of the double summation with respect 
to i and j for a single body of revolution the zero -lift wave drag of 
equation ( 27 ) can be expressed as 

I i-x 

Di * U ^ ^/fi"(x i )Z^"(xj)l 1 (x i ,Xj) (28) 

i=2 j=i 

Multiple-component configurations . - For configurations whose 
equivalent bodies of revolution satisfy the end conditions given by 
equation (4), the zero-lift wave drag of each equivalent body of revolu- 
tion is given by equation (28). The zero-lift wave drag of the given 
configuration is determined by use of the results of equation (28) for 
each equivalent body of revolution in equation (l). If the end conditions 
of equation (k) are not satisfied, the results of equation ( 15 ) for each 
equivalent body of revolution (obtained from actual oblique areas) must 
be used in equation (l). It should be noted that equation (24) cannot 
be applied to calculate the total zero-lift wave drag of multiple- 
component configurations, since this equation is valid only for isolated 
bodies of revolution for which S T (Z) ^ 0. 
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Sheared wing 
for 0 ^ 9 < -j- 


An analysis of the integral equation for a single equivalent body 
of revolution (eq. ( 3 )) reveals that the total zero-lift wave drag of an 
equivalent body of revolution of a multiple -component configuration can 
be calculated as the sum of the zero -lift wave drag of each component 
alone plus the zero-lift interference wave drag between all possible pairs 
of components. For simplicity consider the complete wing shown by the 
solid lines in sketch (b). Each 
wing panel is considered to be one 
component of a multiple -component 
configuration. The dashed lines 
in sketch (b) represent the 
,T sheared M configuration ( des cribed 
in ref. 7 for a given 0 and Q ) 
whose normal area distribution is 
that of the equivalent body of 
revolution. From sketch (b) the 
total second derivative of the 
combined area distribution of the 
two wing panels (i.e., the equiva- 
lent body of revolution) is 

Sketch (b) 


Basic wing and 
sheared wing for 

e-- 



Siji"(x) - Sg 11 (x)+Sg (x) 


and the total zero-lift wave drag of the equivalent body of revolution 
is given by 


/ x a x a 

f [S E "(x 1 )+S H "(x 1 ) ] [S E "(x 2 )+S H "(x 2 )]Zn|x 1 -x 2 | dx!dx 2 

- x b ' X b ( 29 ) 


Expanding the integrand of equation (29) yields 


Dp 


JL 

2tt ( 


a, 

Sg" (x]_ )Sg" (x 2 ) Zn Xi-x 2 dxxdx^- 


-x> 


JL 

2 * 




(x 1 )S H ,l (x 2 )Zn |x x -x 2 |dx 1 dx 2 - 


~ A b 


q 

Tt 



px a 

/ S E ”(x 1 )S H ,, (x 2 )Zn |x!-x 2 JdXidx^ 

-x-u 


(30) 


The first two terms of equation (30) represent the zero -lift wave drag 
alone of each component (wing panels E and H, respectively) of the equiv- 
alent body of revolution of the complete wing. The last term represents 
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the mutual zero-lift interference wave drag between the components of the 
equivalent body of revolution. Equations ( 29 ) and (30) indicate that the 
analysis is performed over the total length of the sheared configuration. 
Since the second derivative of the area distribution of each sheared wing 
panel (sketch (b)) is zero over some portion of the total length, 
equation (30) can obviously be written 


Dip 



S E "( x i)S E n (x 2 )Zn X!-x 2 1 dxxdx^- 


JL 

2rr 



Sjj" (x 1 )Sjj T1 (x 2 )Zn | x x -x 2 1 dx!dx 2 - 


— r Xa 

2 / S E "(x 1 )S H "(x 2 )ln|x 1 -x 2 |dx 1 dx 2 ( 31 ) 

J G ° 

“ x b “ 

Equation ( 31 ) indicates that the analysis is performed over the individual 
lengths of each component of the sheared configuration. The analytical 
solutions of equations (30) and (31) yield identical results; however, 
as will be subsequently demonstrated, the approximate solutions of these 
two equations may differ. 



For convenience, the first term of equations (30) and (31) will be 
designated as Dgg, the second term as Djqj, and the last term as Deh- 
In the notation of the finite series form developed above, the approximate 
expressions for Dgg, an< l -^EH are * for an analysis performed over 

the individual length of each equivalent -body component (eq. (31)) 


!e i-i 

^EE ~ ^ ^ AS E 1, (x Ei )^SE ,! ( x Ej ) I i( x E i > x Ej ) (32) 

i=2 

I H i-i 

%H * 2 ^ ^ ^^H”( x Hj_)^H M ( x Hj ) I i( x H i ^ x Hj) (33) 

i=2 j=i 

!e t H 

DeH ~ 2* Zj ^^E"( x E i ) AS H"( x Hj) I i( x E 1 ^ x Hj) (34) 

1=1 j=l 


For an analysis performed over the total length of the equivalent body 
of revolution, equations (32) to ( 34 ) are applicable when the subscript 
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of the upper limit of summation is dropped and I represents the number 
of points used to divide the total length into equal intervals. Finally, 
the total zero-lift wave drag of the given configuration is, in accordance 
with equation (l). 



Extending this procedure to configurations with more than two components 
the total zero -lift wave drag can obviously be calculated as the sum of 
the zero-lift wave drag of each component alone plus the mutual zero-lift 
interference wave drag between all possible pairs of components. On the 
other hand, if each equivalent body of revolution is treated as a single 
entity and its total zero-lift wave drag is calculated by equation (29 ) 9 
the total zero-lift wave drag of the given configuration is written in 
the form of equation (l) as 


r\ 27T 

Dt( 3) = ^ J »r(P,e)<ie (36) 


Generally speaking, for a given (3, the averaging process of 
equations (35) and (36) is performed graphically because the area distri- 
butions of the components of the equivalent bodies of revolution change 
with G. The term "averaging process" implies integration with respect 
to G and division of this result by 2jt. In the special case of systems 
of slender bodies of revolution, whose individual -component area distri- 
butions do not change with 0, it is expedient to calculate the total 
zero-lift wave drag by equation (35) for an individual -length analysis. 

For this case the zero-lift wave drag of each component alone is calculated 
directly without recourse to averaging since the area distributions of 
all equivalent bodies of revolution for each component alone are identical. 
The interference wave drag between pairs of bodies, on the other hand, 
must be averaged, as the total area distributions of the equivalent bodies 
of revolution for pairs of components do change with 0. As a result of 
the present method for calculating zero-lift wave drag, it is possible, 
in this special case, to perform the averaging process analytically in 
closed form. 


Interference between a pair of bodies of revolution, S f ( -Zi)=S* ( Z )= 0 . - 
Cons ider the general arrangement of a pair of bodies of revolution, E 
and H, whose centers are longitudinally separated by a distance 5 and 
whose axes are laterally separated by a distance d as shown in figure 2. 
Let each of the bodies satisfy the condition of equation (4), and let each 
body be sufficiently slender that the frontal projection of the area dis- 
tribution formed from each body by the oblique planes can be considered 
to be the same as the normal cross-sectional area distribution of each 
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body (i.e., independent of (3). The zero-lift interference wave drag 
between bodies E and H is given by the last term of equation (35). 
Preparatory to developing this term in finite series form, consider the 
area distribution of the equivalent bodies of revolution for the 
configuration of figure 2 . 

The area distribution (and all derivatives thereof) of each equivalent 
body of revolution of the configuration shown in figure 2 can be found in 
terms of sheared configurations in a manner similar to that employed in 
reference 7 ; that is, the area distribution of each equivalent body of 
revolution can be determined as the normal cross-sectional area distribu- 
tion of a properly sheared configuration. For a specified (3 and 0, 
consider the trace in the xy plane of the oblique plane which passes 
through the center of body H (see fig. 2). Using the coordinate system 
of body E, the area distribution of the equivalent body of revolution at 
x E would consist of the frontal projection of the area of body E at 

x-g^ plus the frontal projection of the area of body H at xjj = 0. Within 
the slenderness requirements of the theory, this projected area is the 
same as the normal cross-sectional area of a configuration consisting of 
body E and of body H shown by the dashed lines in figure 2. In other 
words, body H is sheared a distance | (3d cos 0| rearward for 0 < 0 < n/2 
and forward for n/2 < 0 < ir. The total area distribution of the equiva- 
lent body of revolution consists of the sum of the normal cross-sectional 
area distribution of body E plus that of body H properly translated. 
Since the effect of 0 on the area distribution of the equivalent bodies 
of revolution is merely the translation or shearing of one of the original 
bodies, the second derivative of the area distribution of body H can be 
determined relative to the coordinate system of body H and related to 
the coordinate system of body E by the transformation (see fig. 2) 

X H(E) = x H +5+ P d cos 9 (37) 

When equations ( 3 I+) and ( 37 ) are combined, the last term of 
equation ( 35 ) becomes 

f I E % 

^AS E "(x Ei )^6H"( x Hj )li[ x E i ^ x H(E)j l d0 1 

i =1 

(38) 

where from equations ( 17 ) and ( 37 ) 




2tt 


2n , 


U( x Ej > X H(E) j ) 


= I i[ x E 1 ^ x H(E) J ^ d ^J 

= (x E . -Xjj.— 8— |3d cos 0) 2 Zn x E -Xu -5-(3d cos 0 
1 J I J 


To simplify the notation, equation (39) is written 


( 39 ) 
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I 2 (aj_j ,3d,0) = (a^j-pd cos 0) 2 7n |ajj-f3d cos Q | (40) 

where 

a ij = l x E i -( x H J + 5 ) I (^i) 

For both bodies AS"(x) is independent of Q and it can be determined 
with the aid of figure 2 that the sheared configurations are symmetrical 
for 0 < 6 < it and it < 6 < 2it; therefore equation (38) can be written 

X E % 

Deh( 34) - ^ £ ^AS E "(x E .)^ H "(x H .)l(a i j,pd) (42) 

i=i j=i 


where 


Ka^^d) = i 


n 

I2(aij,pd,0)d0 


o 


(U) 


A detailed evaluation of equation (43) is given in closed form in 
appendix B in terms of the dimensionless coordinates described in 
appendix A. 


Equation (42) is also applicable to nonlifting slender bodies of 
revolution whose equivalent -body area distributions are not considered 
identical to the normal cross-sectional area distribution of the body, 
but are dependent upon (3. For these cases l(aij,3d) of equation (42) 

is unchanged but AS"(x) must be replaced by AS"(x,3) where AS"(x,3) is 
obtained for each body from the second derivative of the frontal projection 
of the area distribution formed from each body by the oblique planes. 


The importance of equation (42) is that for each 3d it yields the 
zero -lift interference wave drag between a pair of bodies of revolution 
from a single set of calculations by an electronic computing machine 
(punch cards). Thus, with tabulated values of equation (43), a single 
set of calculations evaluates each term of equation (35)* On the other 
hand, to evaluate equation ( 36 ) a separate set of calculations for each 
value of 6 and an integration are required. Consequently, it is imme- 
diately apparent that the total zero-lift wave drag of a system of slender 
bodies of revolutions can he calculated more rapidly and more accurately 
if the system is subdivided into components and an individual -length 
analysis is used. 

Equation (42) is valid for all 3* &} and 5 greater than or equal 
to zero. Several interesting results are indicated by equation (42) when 
combinations of these variables are zero. If 3 and/or d equals zero 
the averaging process is automatically eliminated (see eqs. (40) and (43)). 
In this case equation (42) yields the same result as equation (34). In 
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particular, if pd = 0 by virtue of the fact that d = 0 the body axes 
are coincident and for this condition equation (42) is applicable to each 
pair of components of the equivalent bodies of revolution of airplane-type 
configurations. For the very special case of two identical bodies of 
revolution arranged so that d = 5 = 0, the result of equation (42), for 
any 3 > 0, is exactly twice that of equation (27); that is, the zero- 
lift interference wave drag, at any supersonic Mach number, between a 
pair of identical bodies of revolution whose axes are superimposed is 
twice the zero -lift wave drag of either body alone. A bar over the sub- 
scripts denote pairs of bodies of revolution which satisfy these 
conditions; thus 

Dgjj(p) = 2Dgg(p) = 2Dgjf(p) (44) 

Equation (44) is also valid for a pair of identical bodies of revolution 
which satisfy the conditions p = 6 = 0 and d f 0, and for the equivalent 
bodies of revolution of a complete wing for p > 0 and 0 = it/ 2 (see 
sketch (b)). The latter result is immediately obvious from equation (3l)> 
since, for 0 = jt/2, S^ ! ’(x) = Sjj t, (x), -x^ = -c 0 /2, and x a = x G . 

Interference between a pair of bodies of revolution, 

S f (-Zi) =0, S'(Q 4 0 .- It may seem logical that the same procedure 
employed above for a pair of bodies of revolution with S T ( -Z i) =S* ( Z ) =0 
can be used to develop a zero -lift interference wave -drag equation for a 
pair of bodies of revolution with S r (-Zi) = 0 and S ] (z) ^ 0 merely by 
using equation (37) to evaluate j(xj) of equation (24) in analytical 
form for all 0 < 0 < jt. Unfortunately, this is not the case. It will 
be recalled that equation (24) was developed from slender-body theory 
for an isolated body of revolution. However, a zero-lift interference 
wave -drag equation for bodies of revolution with these end conditions 
can be developed by using equation (37) an d the proper limits of integra- 
tion for each body in equation (7). Following the same analysis as that 
used for equation (7); an equation similar to equation (15) will evolve. 
The new expressions for K x , K 2 , and K3 will be a function of Zg, Zjj, 

Xj_, x j , pd, 5, and cos 0. The last term will be identical to the unaver- 
aged expression of equation (38). Each term of the new equation can be 
averaged analytically in closed form. The averaged expression for the 
last term of the new equation will, of course, be given by equation (42). 
It is to be emphasized that for bodies of revolution with S'(-Zi) = 0 
and S f ( Z ) ^ 0 it is mandatory to use the second derivative of the fron- 
tal projection of the area distribution formed from each body by the 
oblique planes . 


Procedure for Applying the Method 


Zero-lift wave drag .- Application of the equations developed above 
depends upon the geometry of the configuration being considered. The 
equations of this report to be used to calculate the zero-lift wave drag 
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of each component of an equivalent body of revolution and the zero-lift 
interference wave drag between pairs of components of an equivalent body 
are summarized in table I for several configurations together with the 
averaging process generally required to obtain the zero -lift wave drag 
of the given configuration. 

As a result of the basic approximation of S"(x) each of uhe 
equations listed in table I is easily adapted to electronic -machine com- 
puting techniques, in particular, to punch-card computing machines. 

Values of I 1 (xj_,xj) of equation (28) and of l 2 ( x Ej_; x Hj) equation (34) 

can be obtained from values of l(a^j,|3d) of equation (42) for [3d = 0. 
Hence, if the present method is to be adopted, cards for the values of 
l( a i j , Pd) can be punched for many values of the arguments and permanently 
stored for reuse with values of AS"(x) of a particular configuration to 
be analyzed. This accomplished, the problem of the application of the 
present method resolves itself into the question of how the values of 
2£> T 1 (x) at the points which divide S"(x) into equal intervals are to be 
recorded on punch cards. Generally it will not be possible, or even 
feasible, to calculate AS n (x) at each of these points. Therefore, the 
procedure employed is to record the values of S"(x) on a set of punch 
cards with the aid of an electronic machine called a telereader ; that is, 
one can determine S"(x) at the desired points from a plot of S"(x) vs. x. 
Then, as noted earlier, 


AS’^Xi) = S n (x i )-S"(x i _ 1 ) 


This procedure has been found to yield numerical solutions of the 
zero-lift wave-drag equations in very good agreement with analytical 
solutions of equation ( 3 ) for configurations with no singular points in 
S"(x). When singularities occur in S"(x), the telereader yields values 
of S"(x) in the region of the singular points which are subject to 
errors. These errors have been found to have a predominant influence on 
the calculated wave -drag results when singularities occur at the ends of 
the distribution. When singular points in S"(x) occur in the inboard 
region of the distribution these errors have been found to have only a 
small effect on the accuracy of the calculations. In order to circumvent 
the detrimental effect of singularities at the ends of S" (x), it was 
reasoned that, since S"(x) was approximated by constant values over equal 
intervals, S'(x) is linearly approximated over these intervals and could 
be determined by the telereader from plots of a continuous variation of 
S' (x) vs. x. Then it becomes a simple matter to define a mean value of 
the second derivative of the area distribution as 


S m "(x) = 


S'(x i+ i)-S'(xi) 

x i+i -x i 


for Xj_ < x < x 


■1+1 


(^ 5 ) 


from which Z£ m "( x ) is easily obtained. The results of many calculations, 
not presented here in the interest of brevity, demonstrate that the use 
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of S m " ( x ) rather than S"(x) yields zero-lift wave-drag results which 
are invariably more accurate for distributions with singularities. The 
use of S m "(x) rather than S"(x) for distributions with no singularities 
reduces the accuracy of the calculations by less than one tenth of a 
percent. Consequently, it is recommended that S m "(x) always be used in 
the present method. All calculations of this report made by the present 
method have employed ^"(x). 

A final point of procedure in applying the present method concerns 
the arbitrary choice of the number of points selected to divide the second 
derivative of the area distribution into equal intervals. Obviously, the 
number of points selected will affect the accuracy of the zero-lift wave- 
drag calculations in a manner similar to that which the choice of the 
number of terms of the infinite series affects the accuracy of the results 
obtained by the method of references 2 or 3 . It has been found in the 
present analysis and in reference 3 that 201 points (200 equal intervals) 
is a practical limit for which punch-card computing machines should be 
programed. The results of many calculations indicate that in some special 
cases as few as 101 points, or even less, yield zero-lift wave-drag results 
with acceptable accuracy. However, the most accurate results have been 
found to be consistently obtained with 201 points. Consequently, all 
calculations used to evaluate the present method have employed 201 points . 
The results of a special detailed analysis of the effect of the choice of 
the number and the location of equally spaced points on the zero -lift wave 
drag of one configuration are presented in the next section of this report. 

The development of the present method for calculating zero-lift wave 
drag and the outline of the procedure for applying the method have merely 
presumed that the second derivative of the area distribution (S M (x) or 
Sm"(x)) is available for all components of all equivalent bodies of revo- 
lution of a configuration. The actual determination of these distributions 
is probably the greatest task required in order to perform zero-lift wave- 
drag calculations by the present method. An analytical procedure for 
computing S(x) for wings and tail surfaces is available in reference 7 
and has been expanded in reference 8 for wings (with arbitrary thickness 
distributions) in the presence of bodies. A method for finding S(x) 
and S'(x) (from which ^"(x) is obtained) is given in reference 4 for 
bodies of revolution and wing and tail -surface components of equivalent 
bodies. Because the area distribution of wing and tail-surface components 
of equivalent bodies of revolution changes with (3 and 9 , and because the 
airfoil-section thickness distribution is not generally available in 
analytical form, the graphical method of reference 4 can become tedious. 

The need for a more simplified method for finding S(x), S'(x), and even 
S"(x) for wings and tail surfaces, obviously, exists. To satisfy this 
need a simplified method has been developed for the analytical calculation 
of S(x), S' (x) , and S"(x). The formulation of analytical expressions 
for these distributions obviously depends upon an analytical expression 
for the plan-form thickness distribution. Analytical approximations of 
the airfoil-section thickness distribution as well as the complete 
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plan-form thickness distribution are given in appendix C. The analytical 
expressions are given in terms of a single parameter which identifies 
each of the sheared configuration plan forms described in references 7 
and 8. This single parameter combines the dependence of the sheared plan 
forms upon both 3 and 9. The analytical method for calculating S(x), 
S’(x), and S"(x) is presented in appendix D. 

Wave drag.- The general form of equation (15) can be employed to 
calculate the linear-theory wave drag for planar and nonplanar lifting 
configurations when, in addition to the second derivative of the area 
distribution, the force term of equation (2) is known. When the force 
term is also known the second derivative of the area distribution of 
each equivalent body of revolution is represented by the combined term 

S"(x,p,0) - L'(x,p,0) 


and for a specified 3 and 9 equation (2) has ’jhe same form as 
equation (3). Hence in equation (15) one merely replaces AS"(x) by 


A 


S"(x) 


-s L ' (x) 


If it is desirable to separate the effects of thickness and lift on the 
wave drag, equation (2) can be expanded by multiplying out the integrand. 
This manipulation yields three terms for the wave drag of each equivalent 
body. Two of the terms, identical or similar in form to equation (3), 
represent the separate contributions to the wave drag of thickness and 
lift and can be evaluated using the form of equation (15). The third 
term, similar in form to the last term of equations (30) and ( 31 ), repre- 
sents the interference wave drag between thickness and lift and can be 
evaluated using the form of equation (42) with values of l(afj,pd) for 
3d = 0. In general, L'(x,3,0) is difficult to obtain; however, a review 
of the literature revealed that reference 9 provides a systematic method 
for computing wing pressure distributions for wings with subsonic leading 
edges and supersonic trailing edges. From this information it is possible 
to obtain L(x,3,0) and L' (x,P,0) . A procedure for finding these quan- 
tities in terms of the lift distribution of sheared wings is briefly 
outlined in appendix E for wings with zero taper ratio. 

Vortex drag.— The general form of equation (15) can also be used to 
calculate the vortex drag of a finite -span wing. The familiar form of the 
equation for the vortex drag is 

b b 

% = - A f * f* r’(y 1 )r , (y 2 )ln|y 1 -y 2 idy 1 dy 2 (46) 

~ 2 ~ 2 
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where p is the free-stream density and r r (y) is the slope of the 
spanwise circulation distribution. Equation (46) is identical in form 
to equation (7). By means of the same approximation of T T (y) as was 
used for S"(x), the vortex drag can be calculated by equation (15) with 
q, AS" (x) , x, and Z replaced by p/2, A T T (y), y, and b/2, respectively. 


EVALUATION OF THE METHOD 


The method will be evaluated by calculating the zero-lift wave drag 
of several bodies of revolution and comparing the results with those 
obtained by analytical solutions (when available) and by the method of 
reference 3. The calculations are restricted to Mach number 1 in order 
to avoid the integration required of the method of reference 3 at Mach 
numbers greater than 1. Subsequently, the merits of the present method 
for Mach numbers greater than 1 will be briefly examined. The two basic 
arrangements of the bodies considered in the evaluation are presented in 
figure 3- As shown in figures 3( a ) and 3(B) the body centers are separated 
longitudinally a distance 5 = 0 (8 = 0) and 6 = 2Z (6=2), respectively. 
For these arrangements the wave-drag analysis is performed over the total 
lengths 2Z and 4Z, respectively. For the arrangement of bodies with 
6=0 the wave -drag analysis is also performed over the alternate length 
4z as indicated in figure 3(c). Equations defining the dimensionless 
area distribution of the bodies used in these arrangements are given below: 


Body 

Equation of area distribution 


1 

3i(s) = | a -ef >2 

(47) 

2 

S 2 (i) = | a -I 2 ) 572 

(48) 

3 

s 3 (i) = || U-i 2 ) 7/2 

(49) 

1 + 

s 4 (i) =V? d-^) 9/2 
63 

(50) 

5 

§ 5 ( 1 ) = | id -i 2 ) 3/2 

(51) 

6 

s 6 (i) = | id-i 2 ) 5/2 

(52) 

7 

Sr(l) = § i(l-i 2 f /2 

(53) 

8 

s 8 (i) = ^ i(i-i 2 ) 9/2 

(54) 


For each body -1 < g < 1. Equations (1+7) to (50) are even functions with 
the coefficients so chosen that the dimensionless value of the volume of 
each body is it/2. Equations (51) to (54) are odd functions corresponding 
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to bodies 1 to 4, respectively, and the volume of each of the bodies 5 
to 8 is obviously zero. It should be noted that bodies 1 and 5 have 
singularities in §"((■) at both ends (i.e., at £ = -1 and f- = 1). 

Comparisons will be made of the total zero-lift wave drag for 
several pairs of bodies, the separate contribution of each body alone to 
the total, and the separate contribution of the interference between 
bodies. The interference drag results will be considered first, since, 
as discussed in the METHOD section and amplified in the following para- 
graphs, these results provide information from which calculations can be 
made of both the total wave drag of a pair of bodies and the contribution 
of each body to the total. Before presenting the results of the calcula- 
tions it is desirable to mention the specific computing procedures 
employed for the bodies considered herein. These specific computing pro- 
cedures were employed in order to obtain consistently the most accurate 
results by each method. 

Only the zero -lift interference wave drag was calculated by an 
electronic computing machine for both the present method and the method 
of reference 3. Values of 3(0 and S m n (|) were calculated at each of 
the 201 points from equations (47) to (5*0. The existing machine - 
computing program at Ames Research Center for the method of reference 3 
yields directly, for a given value of 3 cos G and up to 49 terms, only 
the total zero-lift wave -drag results based on an analysis performed 
over the total length of a configuration. However, interference drag 
results were obtained by a special manipulation of the data obtained from 
the existing program. This manipulation is described in appendix F which 
presents the computing procedure used in this report for the method of 
reference 3 . 

Four significant figures were retained for all machine -computed 
calculations of the interference drag results for both methods. For the 
configurations of figure 3(*b) with singularities in S"(g) at the ends 
of at least one body, only two significant figures were possible In the 
case of the present method. This loss of significant figures is a direct 
result of the large values of S m M (0 at and near the ends of the body. 

The total zero-lift wave drag of a pair of bodies of revolution is 
given by the sum of the contribution to the total of each body alone and 
the interference between bodies (see eq. (35)). The body-alone contri- 
bution can be determined from the interference results in accordance with 
equation (^4). Hence, by combining equations (35) and (44) the following 
equations can be written for the dimensionless total zero-lift wave drag: 

- jj Sg^(a) -i- ^ Dj^(a )+Dgjj(a ) (55) 


for the configurations of figure 3( a )j 
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- \ + \ Drtt( c ) + %Tt (^) (56) 

for a total-length analysis of the bodies in figure 3(b), and 

%■ = \ + \ ( 57 ) 

for an individual-length analysis of the bodies in figure 3 (b). In 
equations (55) to (57) the terms in parentheses designate that part of 
figure 3 which indicates the configuration length over which the wave- 
drag analysis is performed. 


Accuracy of Calculations for Interference 


The dimensionless zero-lift interference wave drag for the bodies 
in the arrangements of figure 3 is presented in table II. Rough sketches 
of figures 3(a), 3(b), and 3(c) are duplicated in the column headings to 
indicate the arrangement of the bodies and the length over which the wave- 
drag analysis is performed. In column 1 of table II are tabulated all 
possible combinations of pairs of bodies whose area distributions are 
defined by equations ( 47 ) to ( 54 ). Note that the bodies with singularities 
in S"(£) at the ends are marked with an asterisk. Columns 2 and 3 of 
table II list the analytical values of the wave drag. For 8=2 the 
analytical solution of the last term of equation (31) for nonidentical 
bodies has not, as yet, been obtained; therefore, analytical solutions 
for only identical pairs of bodies are presented in column 3. However, 
it is felt that the pairs of identical bodies with 8=2 are sufficient 
to evaluate the accuracy of the present method for calculating zero-lift 
wave drag. 

Bodies with 8 = 0 . - For all pairs of bodies in the arrangement of 
figure 3(a) for which data are available both methods yield zero -lift 
interference wave-drag solutions which are accurate to within a fraction 
of 1 percent of the analytical solutions (compare columns 2, 4, and 7 of 
table II). For configurations with no singularities or discontinuities 
in §"(|) the method of reference 3 yields exactly the analytical solu- 
tions. The data of columns 5 and 8 of table II are included merely as 
data from which the body-alone contribution to the total zero-lift wave 
drag can be obtained for configurations of figure 3(b) which are analyzed 
over the total length of the configurations. It will be noted that ana- 
lytical solutions for configurations of figures 3(a) and 3(c) are identical 
(see column 2 of table II); whereas, the numerical solutions by both 
methods are less accurate for the configurations of figure 3(c). 

Bodies with 5 = 2 . - For all pairs of bodies in the arrangement of 
figure 3(b) for which data are available both methods yield zerc-lift 
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interference wave -drag solutions which are accurate to within a fraction 
of 1 percent of the analytical solutions for pairs of bodies with no 
singularities or discontinuities in S"(|) (compare columns 3, 6, and 9 
of table II). For pairs of bodies, at least one of which has singularities 
in S n (£) at the ends, the method of reference 3 and the present method 
yield zero-lift interference wave-drag solutions which are accurate to 
within -8 and +14 percent of the analytical solutions, respectively. 
However, these figures are percentages of small values of the interference. 
Furthermore, it should be recalled that with the use of the existing 
machine -computing program the results for these bodies by the present 
method are based on only two significant figures. Bodies 1 and 5, for 
which calculations by the present method are in greatest error, are bodies 
with singularities at the ends which have minimum zero-lift wave drag at 
transonic speeds for a given volume and length and a given first moment 
of area, respectively. For these reasons an evaluation of the accuracy 
of the zero -lift wave -drag results computed by the present method for 
these bodies should constitute the most severe test of the method. 


Accuracy of Calculations for Body Alone 


In view of the simple relationship given by equation (44) between 
the zero-lift wave drag for a body alone and the interference between a 
pair of identical bodies with 6=0, body-alone results have not been 
tabulated. It is obvious, however, that both methods yield numerical 
solutions for the body-alone zero-lift wave drag which are accurate to 
within a fraction of 1 percent of the analytical solutions. 


Accuracy of Calculations for Total Configurations 


The dimensionless total zero-lift wave-drag results for only pairs 
of identical bodies in both basic arrangements of figure 3 have been cal- 
culated using the data of table II in equations (55) ho (57) and are 
presented in table III, These pairs of bodies are considered sufficient 
to evaluate the accuracy of total zero-lift wave-drag calculations by the 
present method in that these bodies represent the extreme deviations from 
the analytical solutions and are generally typical of all the pairs of 
bodies in their demonstration of the effect of the different lengths used 
in the analysis. The data in parentheses in table III indicate the percent 
deviation from the analytical solutions. 

Bodies with 6=0 .- Since the total zero-lift wave drag for these 
bodies was obtained from the zero -lift interference wave drag of these 
bodies, it is clear that the accuracy of the total wave-drag results for 
these bodies is consistent with that of the interference wave -drag results; 
that is, for all pairs of bodies in table III both methods yield total 
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zero-lift wave-drag solutions which are accurate to within a fraction 
of 1 percent of the analytical solutions (compare columns 2, k, and 5 of 
table III). For configurations with no singularities or discontinuities 
in S ,T ( £ ) at the ends, the method of reference 3 yields exactly the 
analytical solutions. 

Bodies with 5=2 . - Total zero -lift wave drag for these bodies can 
be calculated from an individual- or a total-length analysis. As noted 
in the METHOD section and demonstrated in column 3 of table III, the 
configuration length used in the analysis does not influence the analyti- 
cal solutions. The length does, however, affect the numerical solutions. 
The effect is small for the simple bodies of revolution of this report, 
even for bodies with singularities at the ends of S"(|) (see table III). 
Preliminary calculations for a more general airplane -type configuration 
indicate, however, that the influence of the length may not be small for 
such configurations in that S n (£) of the sheared wings of many equivalent 
bodies is discontinuous. Therefore, the authors feel that some discussion 
of the length is warranted for the present configurations even though, in 
this case, the effect is small. 

For pairs of bodies with no singularities or discontinuities in 
§"(£) both methods yield results which are accurate to within one third 
of 1 percent of the analytical solutions regardless of the length of 
analysis (compare columns 3* 1> and 9 for the present method and columns 3> 
6, and 8 for the method of ref. 3)- Even though the deviations from the 
analytical solutions are very small for each length of analysis, a defi- 
nite trend toward more accurate solutions for the individual -length anal- 
ysis is evidenced by the results for these bodies. For pairs of bodies 
with singularities in §"(£) at the ends (bodies lxl and 5x5) both methods 
yield results which are accurate to within ±3 percent of the analytical 
solutions. The same trend toward more accurate solutions for the 
individual -length analysis is evidenced by the results for these latter 
bodies for the method of reference 3 but not, however, for the present 
method. In spite of these seemingly anomalous results for the present 
method, and even though the effect of length of analysis is small, the 
data of table III indicate that both methods yield comparable total zero- 
lift wave-drag solutions which are in closer agreement with analytical 
solutions when the wave-drag analysis is performed over the individual 
lengths of the bodies. 

Bodies with 0 < 5 < 2 and 6 > 2 . - The total zero-lift wave -drag 
calculations at Mach number 1 (or at any supersonic Mach number) for 
pairs of bodies arranged with these values of 8 will be more accurate 
than those discussed above for 8 = 2 if the wave-drag analysis is always 
based on an individual- length analysis. The body-alone contribution to 
the total zero -lift wave drag is, of course, independent of 5; but both 
the accuracy and magnitude of the interference contribution are affected 
by S. As 8 becomes less than 2, the interference calculations become 
more accurate as a result of a larger number of points of the 201-point 
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analysis lying on each body. As the value of § increases beyond 2, the 
interference calculations are subject to greater inaccuracies , yet the 
magnitude of this portion of the wave drag diminishes to an almost 
negligibly small percentage of the total zero -lift wave drag. 


Accuracy of Calculations at Mach Numbers Greater 
Than 1 for Total Configurations 


A complete evaluation of the present method as a tool for calculating 
total zero-lift wave drag (eq. (35)) should include an evaluation of the 
accuracy of the method at Mach numbers greater than 1 for the configura- 
tions of this report as well as the more general airplane-type configura- 
tions. No quantitative evaluation of this aspect of the accuracy will be 
made; however, from the information available in this and the previous 
section of this report it is possible to make a qualitative evaluation. 

Pairs of bodies of revolution. - The area distribution of each body 
does not change with 0; therefore, the zero-lift wave drag of each body 
alone is obtained, with comparable accuracy, from a single set of calcu- 
lations by each method. The area distribution of the combined bodies of 
each equivalent body of revolution dees change with 0 . For the method 
of reference 3 'the zero-lift interference wave drag between bodies must 
be computed for several values of G and the results averaged as dis- 
cussed relative to equations (35) and ( 36 ); whereas, for the present 
method the zero -lift interference wave drag is obtained from a single set 
of calculations, analytically averaged for all 6 . As a result of the 
averaging process required of the method of reference 3> it can be stated 
that the present method yields the mere direct and probably the more 
accurate total zero-lift wave-drag results for systems of bodies of 
revolution. 

Airplane-type configurations.- The area distribution of the various 
components of airplane -type configurations changes with 0. Hence, for 
these configurations both methods, in general, yield total zero-lift wave- 
drag results with a comparable accuracy. If a particular configuration 
contains several body-of -revolution components such as the airplane-type 
configuration shown in table I, the present method will probably yield 
more accurate results for the reasons mentioned above concerning the 
zero -lift interference wave drag between bodies of revolution. 


Computational Time Required 


An evaluation of the computational time required by the present 
method is equally as important to the complete evaluation of the method 
as is an evaluation of the accuracy of the method. Just as the most 
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accurate zero-lift wave-drag values obtainable by the present method and 
the method of reference 3 were used to effect an equitable evaluation of 
the relative accuracy of the methods, so it is desired to evaluate the 
relative computational time for both methods on the basis of equitable 
electronic-computing machine techniques. Accordingly, this evaluation 
will be based upon the following assumptions: 

i. The zero-lift wave-drag equations for an equivalent body of 
revolution for both methods are programed in such a manner to permit 
calculations of the total zero -lift wave drag of multiple -component con- 
figurations from separate calculations of the zero-lift wave drag of each 
component alone plus the zero-lift interference wave drag between 
components . 

ii. The initial data for each method are obtained by a telereader 
from a plot of S(|) or S'(|) vs. £, as the case may be. It should be 
recalled that values of S m "(|) are calculated from values of S'(£) 
according to equation (53).> and this operation is assumed to be included 
in the program. 

iii. S(|) and S'(|) are obtained for wings and tail surfaces by an 
analytical method similar to that presented in appendix D. 

The present method and the method of reference 3 require the same 
computing time to calculate the total zero-lift wave drag of any configu- 
ration at Mach number 1. If plots of the distributions are available, 
calculations for one equivalent body of revolution require up to 1-1/2 
hours from the beginning of the telereader operation to the final answer. 
For configurations with components whose equivalent body area distributions 
change with 9 , both methods require the same computational time because 
the time required to calculate the total zero-lift wave drag of each 
equivalent body is identical and the averaging process can be performed 
similarly for both methods. For configurations consisting entirely of 
body-of -revolution components, both methods require the same time to cal- 
culate the zero-lift wave drag of each body alone. The present method, 
however, affords a saving in time required to calculate the zero-lift 
interference wave drag between pairs of bodies because the averaging 
process cannot be performed similarly for both methods. The present 
method inherently averages the interference contribution of all pairs of 
equivalent bodies with a single set of calculations; whereas, by the 
method of reference 3; "the averaging process must be performed for the 
results of a finite number of pairs of equivalent bodies of revolution. 
Consequently, the saving in time is directly proportional to the number 
of pairs of equivalent bodies (or values of 0) necessary to define the 
integrand of the last term of equation (35). For example, for two bodies 
of revolution (the minimum required for interference) a minimum of five 
equivalent bodies is required. Hence, zero-lift interference wave-drag 
calculations by the present method would only require 20 percent of the 
computing time required of the method of reference 3, and the total 
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zero-lift wave-drag calculations would require 60 percent of the computing 
time required of the method of reference 3- If more than five equivalent 
bodies were required, the saving in computing time would be even greater, 
as would be the case for a system of more than two bodies of revolution. 


Ultimate Capabilities of the Present Method 


The analyses presented in this report have developed a numerical 
method for calculating zero-lift wave drag which employs a knowledge of 
the second derivative of the area distribution of a configuration and 
which is easily adaptable to electronic-machine computing equipment. The 
ultimate desire, however, is the development of a method for computing 
zero-lift wave drag by ordinary desk computing techniques. As a result 
of a detailed analysis of the effect of the choice of the number and the 
location of the points used to divide S m "(£) into equal intervals, such 
a method appears possible for components of a configuration. Presented 
in figure 4 are the results of this detailed analysis for a single body 
of revolution whose area distribution is composed of three parabolic arcs. 
The dimensionless zero-lift wave drag tabulated in figure 4 demonstrates 
that as few as 25 points yield a solution accurate to within 1.75 percent 
of the analytical solution. This suggests that by a judicious choice of 
intervals (not necessarily equal) a fewer number of points may be con- 
sidered which will yield numerical solutions in almost exact agreement 
with analytical solutions and thereby make desk calculations feasible. 
Preliminary calculations indicate that zero-lift wave-drag values of any 
one component of a configuration which are essentially identical to the 
analytical values can be obtained with as few as 10 unequal intervals if 
the intervals are systematically determined from information concerning 
the source strengths which define the geometry of the component. 


CONCLUDING REMARKS 


A numerical method based on linearized and slender-body theories, 
which is easily adapted to electronic -machine computing equipment, has 
been developed for calculating the wave drag of single- and multiple - 
component configurations from a knowledge of the second derivative of the 
area distribution of a series of equivalent bodies of revolution. Zero- 
lift wave -drag results have been calculated at a Mach number of 1 for 
several simple analytical shapes (bodies of revolution) by the present 
method and by a method which employs the Tchebichef form of harmonic 
analysis of the configuration area distribution. The relative accuracy 
and computational time required of both methods to calculate zero-lift 
wave drag at a Mach number of 1 have been evaluated by comparing the 
zero-lift wave -drag solutions for these simple shapes computed by both 
methods with those obtained from analytical solutions of the wave-drag 
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equations. From these results it was possible to evaluate also the 
relative accuracy and computational time required of both methods to cal- 
culate the zero-lift wave drag at supersonic speeds. The following remarks 
are warranted as a result of these evaluations. 

For the pairs of bodies considered in this report at a Mach number 
of 1, both methods yield numerical solutions for the zero-lift wave drag 
of each body alone to an accuracy of a fraction of 1 percent of the ana- 
lytical solutions. Calculations by both numerical methods of the zero- 
lift interference wave drag between bodies were generally accurate to 
within a fraction of 1 percent; however, in some special cases for which 
the interference was small, these results were accurate to within only 
-8 and +14 percent of the analytical values for the harmonic analysis and 
the present method, respectively. Total zero-lift wave-drag solutions 
generally were calculated more accurately as the sum of the zero -lift 
wave drag of each body alone plus the interference between bodies. All 
total drag results for bodies with no singularities or discontinuities in 
the second derivative of their area distributions were accurate to within 
a fraction of 1 percent of the analytical solutions, and for bodies with 
singularities the results were accurate to within 3 percent of the 
analytical solutions. Both methods required the same computing time. 

For total zero -lift wave -drag calculations at supersonic Mach numbers, 
the relative accuracy and computing time required by both methods is com- 
parable for airplane-type configurations. For systems of bodies of revo- 
lution both methods yield similar results with comparable accuracy; however, 
the present method only requires up to 60 percent of the computing time 
required of the harmonic-analysis method for two bodies of revolution and 
less time for a larger number of bodies. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Oct. 17, 1958 
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APPENDIX A 


DIMENSIONLESS ZERO -LIFT WAVE -DRAG EQUATIONS 


In order to take advantage of the decrease in the number of parameters 
resulting from similarity considerations and to facilitate the calcula- 
tions, the quantities defined in the METHOD section will be made dimen- 
sionless. The procedures presented apply to all of the sheared or 
equivalent configurations previously discussed. However, in the interest 
of simplicity, the dependence of the quantities upon 3 and 0 will be 
omitted from the notation. 


To establish a general dimensionless coordinate system each coordinate 
is divided by a characteristic length, areas are divided by some reference 
area, volumes are divided by some reference volume, and so on. In this 
manner the relationship between dimensional and dimensionless lengths and 
longitudinal area distributions can be written 


I 



n 



i 



S(i) = S(x) 
S R 


for 


-Z i < x < Z 

_ ii < f < _L 

. k Y k Y 


(Al) 

(A2) 

(A3) 

(AM 


where k x , ky, and k z are characteristic lengths in the x, y, and 
z directions, respectively, and Sp is the reference area. From a 
differentiation of equations (Al) and (a 4) the relationship between the 
derivatives of the dimensional and dimensionless area distributions are 
found to be 

S>(|) = jg S' (x) (A5) 

S"(D = hLl s"(x) (a6) 


With the aid of equations (Al) to (A6), the general relationship between 
the dimensional and dimensionless zero -lift wave drag can be shown to be 
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D 



(AT) 


Values of the characteristic lengths and reference areas are tabulated 
below for bodies of revolution and wings or tail surfaces* 


Configuration 

k x 

k y 


Sr 

Body of revolution 
Wing or tail surface 

l! = l 

c 0 /2 

*o/2 

t Q 

So 

^o(^o/2) 


z 



The characteristic length l or h 
is the half -body length and S Q is 
the maximum frontal area of the 
body. The characteristic lengths 
for wings or tail surfaces are 
shown in sketch (c). By means of 
equations (Al) to (A6 ) and the 
above table , the dimensionless zero- 
lift wave -drag equations for both 
the present method and the method 
of reference 3 can determined 
as discussed in the following 
paragraphs . 


PRESENT METHOD 

Single Body of Revolution, S f (-Zi) =0, S 1 ( Z ) / 0 


The zero-lift wave drag for this kind of body is given by equation (5); 
from which the dimensionless form is determined to be 


D = - 


2jt 


i 


3' , (6 1 )3"(! 2 )(ln|S 1 -§ 2 |+Zn z)d^d| 2 + 


-1 -1 


1 [S' (1) 

7T 1 2 


, !n isfiy • lniJT °_ 


S' (1) J S"( | 2 )(ln |l-| 2 |+in i)d| 

"I 

(A8) 
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In accordance with equations ( 7 ) and (8), the dimensionless expressions 
T) 1 and D 2 are given by the first and second terms of equation (A8), 
respectively. 


Single Body of Revolution, S*(-Zi) = S f (z) = 0 


For this kind of body the second term of equation (A8) is identically 
zero and the dimensionless zero-lift wave drag is given by only the first 
term. This first term can be simplified, however, for bodies of revolution 
with S T (-l) = S T ( 1 ) = 0. After an integration by parts, it is a simple 
matter to show that the Zn l term vanishes. Hence, 

Bi = f S ,, (l 1 )§"(l 2 )zn|l 1 -6 2 |di 1 d5 2 (A9) 

-1 -1 


Multiple -Component Configurations 


The general relationship between the dimensional and dimensionless 
equations for an equivalent body of revolution of a multiple -component 
configuration is given by equation (A7). The total zero-lift wave drag 
of an equivalent body of revolution of a two -component configuration can 
be expressed by the unaveraged expression for equation (35). Hence, if 
each component is made dimensionless with respect to its own characteristic 
length and reference area, the relationship between the dimensional and 
dimensionless total zero-lift wave-drag equations is 

1)1 • (Si q6EE + (Si qS ™ * (S) E (Sl qBE » (A10) 

E n hi n 

Experience has demonstrated that in the majority of cases it is more 
convenient to make the complete configuration dimensionless with respect 
to only one of the components. In this case equation (A7) is applicable 
to the total configuration and values of S R and k x are used for only 
one (either) component. 

In view of the foregoing, the dimensionless finite-series form of 
the zero-lift wave-drag equations developed in the METHOD section of this 
report can be obtained if AS"(x), x, and Z are replaced by AS"(£), £, 
and Z/k x , respectively, and if the logarithm terms of the equations for 
D 2 and/or T) x are properly evaluated. As a typical example, the dimen- 
sionless zero-lift wave-drag equation for a single body of revolution 
with S 1 ( -Z 1 ) = S 1 ( Z ) = 0 can be written (see eq. (28)) 



OJ 


4 


I i-i 

81 ' ^ Z Z iS,,(?i)AS " (tj)li<ii ’ ij) 

i=2 j=i 

where, from equation (l?), 

= (t i -6j) 2 in|s 1 -s 2 | 


(All) 


(A12) 


METHOD OF REFERENCE 3 

Use of equations (Al) and (A5) in equations (F4), (F2), and (F3) 
yields the dimensionless zero-lift wave-drag equation for an equivalent 
tody of revolution for the method of reference 3 


where 



(A 13 ) 


(A14) 


and 


I 



cos cp 


(A15) 


SIMILARITY CONSIDERATIONS 


With the definitions of k x , k y , k z , and S R tabulated above, a 
general expression for the wave drag, at any supersonic Mach number, for 
a system of bodies of revolution is 


I>B 



<iI>b P 



(A16) 


where 


®Bpf) 





d0 
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and. for an airplane -type configuration is 


DW = 




(A17) 


where 




ae 


Dividing both sides of equations (Al6) and (A17) by q and a proper 
reference area, one obtains the wave-drag coefficient 


C 


_ 1 _ 

°B S B 




C 




(A18) 

(A19) 


For bodies of revolution S B = S D = jr/Udo 2 , where do is the maximum 
diameter of the body. For airplane -type configurations b /Sy = A, the 
aspect ratio of the wing and t 0 /c 0 = t 0 , the maximum thickness ratio of 
the airfoil section. Therefore equations (A18) and (A19) can be written 


C 


% = 



% 



(A20) 

(A21) 


Equations (A20) and (A21) demonstrate the advantage to be gained from 
similarity considerations which are afforded as a result of wave-drag cal- 
culations made in a dimensionless coordinate system. From equations (A20) 
and (A21), it is clear that dimensionless wave-drag calculations of a 
single configuration, in effect, yield the wave-drag coefficients for an 
entire family of related configurations. The wave -drag coefficients of 
geometrically similar families of bodies of revolution are related by the 
fineness ratio of the bodies, and those for geometrically similar airplane - 
type configurations are related by the aspect ratio and maximum thickness 
ratio of the wings. 
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APPENDIX B 

EVALUATION OF THE INTERFERENCE INTEGRAL FOR 
PAIRS OF BODIES OF REVOLUTION 


In the METHOD section of this report the zero-lift interference wave 
drag between a pair of bodies with S(-Zj.) = S(z) =0, in the arrangement 
shown in figure 2, is given by equation (42). If the complete configura- 
tion is made dimensionless with respect to body E, that is, both x coor- 
dinates are divided by Zg and both area distributions are divided by 
S Eo , equation (42) can be written (see appendix A) 

2 J E X H 

deh(S) “&?(¥) X Z^ E " UE i )A§H " UH j )i(ai j ,g) (B1) 

i=i j=i 


where 


i / (a-y-B cos S) 2 Zn|aij-B cos 0|d0 




B - M 


(B2) 

(B3) 

(B4) 


This appendix will evaluate the integral function defined by 
equations (B2) to (B4). Using the substitution 


7 



(B5) 


one can write equation (B2) as 


i(Sij,B) = 


/ It It 

(y-cos 0) 2 Zn(y-cos 0) 2 d0+Zn B \ / (y-cos 0) 2 d0 
o '"o 

(B6) 


Let 


it 

K(y) = ^ J (y-cos 0) 2 Zn(y-cos 0) 2 d0 


(B7) 
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If the integrand of equation (B 7 ) is grouped as Zn(7~cos 0 ) 2 and 
l/2(7~cos 0) 2 d0, integration by parts yields 


K(r) = (y 2 + i) in I7+1 1 - (y 2 + 




sin^ d 0 + 

cos 0 


2 r i 


sin 2 0 
7-cos G 


d G 


11 f sin 2 

2 7- 


i n 2 0 cos G 


cos G 


d G 


(B8) 


Since 


Sin 2 e _ l-cos 2 0 = 7+cos e _ _r^l_ 


7-cos G 7-cos G 


7-cos G 


and 


sin 2 e cos e = cos , 9 , -co sfe = cos 2 0+7 cos 0+ ( 7 2_i) _ zlzflil 
7-cos G 7-cos 0 7-cos 0 


equation (B8)> after the indicated integration of the simple terms is 
performed^ can be written 


K( 7 ) = ( 7 2 + 


l) !n b + i| + 1(3^ + l) - 0 + l)l/ r-S s~f de ~ 


3 / 2.^ 1 jO 

2 7 y J ft j Q 7-cos 0 


(B 9 ) 


If the integrand of the first integral expression in equation (B 9 ) is 
grouped as 0 and (sin 0/7-cos 0)d0, integration by parts yields 


if 


G sin 0 
7-cos 0 


d 0 = Zn 


r * 

7+1 1 - i* / Zn( 7 -cos 0 


)d0 


(BIO) 


Now 


p ft 

i / Zn(7-cos 0)d0 = Z n ^ (I /I + n/t 2 -!) for I l| ^ 1 


(Bll) 
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The partial derivative with respect to 7 of the left-hand side of 
equation (Bll) is zero for | 7 | < 1 . This demonstrates that the integral 
expression is independent of 7 and hence has a constant value for 
| 7} < 1 . The constant is evaluated as in 2 . Therefore 

Zn| 7+l| +Zn 2 -Zn( | ■y\+ s f/^T) for | y\ > 1 

Zn|7+l|+Zn 2 for j 7 1 <1 

(B 12 ) 

The last integral expression in equation (B 9 ) is evaluated as follows: 




0 sin 6 
7 -cos 6 


d 0 = 


1 

or 


ae 


7-cos 9 


= < 


tan' 


//7 2 - 1 tan i 6 
^®=r \ 7-1 

sJT -7 2 tan i 0-f (1+7 ) 


— iTT 


-IO 

71 


In 


\/l- 7 2 Jl-7 2 tan ~ 0 -(l+ 7 ) 

^ — 1 o 


for J 7 | > 1 


for | 7 < 1 


Therefore 


ae 


s/r 2 -! 


for 7 > 1 


7-cos 0 


= \ tor 7 < - 1 


n/t 2 -! 


0 


for 7 3 s 1 


From equations (B 12 ) and (BI3) equation (B 9 ) reduces to 


(B 13 ) 


[l + |) ‘ (j 2 + ^ ln 2 + (j 2 + Zn( | 71+^72-1) - | 1 7I 


K( 7 ) =< 


I V 372 + + 1 ) Zn 2 


for j 7 I >1 
for 1 7 j < 1 
(B 14 ) 


The second integral of equation (B6) evaluates simply as 


i- J (7- cos 0) S d0 = 72 + i 


(B 15 ) 
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By means 
written 


of equations (B5), (B7), (Bl4), and (B15), equation (B6) is 


I(aij,B) = B qK 


a 


29 

B 


a 


2i 

B 


4 


'n b| 


(bi6) 


where K(aj.j/B) is given by equations (B5) and (Bl4). Equation (Bl6) can 
be further simplified by applying the conditions for bodies of revolution 
with S’ (l) = 0. The terms of the equation resulting from substitution 
of equation (Bl6) in equation (Bl) which are constant or are proportional 
to a-M and a±j 2 will vanish as a result of the conditions imposed for 
all values of a-y by equations ( 25 ) and ( 26 ) for zero area gradient at 
the ends of a body of revolution. Therefore, equation (Bl6) can be written 



(W + |S 2 ) [m(|a i j| + 7a ij 2 -B2)-Zn B] - | | §i j | 


2 -b 2 

l(Sij,B) = < 

for | 

l a ijl 

> ® 


0 f ° r 

g ij 1 

< B 

✓ 




(BIT) 


From figure 2 it can be seen that B is the dimensionless longitudinal 
distance defined by the intersection of the Mach wave (0d cos 9 = 0d 
for 0=0°) with the axes of both bodies. Also, Spj is the dimension- 
less longitudinal distance between pairs of points, one on each_body, used 
in the double summation of equation (Bl). Consequently, l(aij,B) = 0 
for | j | < B has the important physical interpretation that any parts 
of a configuration which lie on or outside the Mach cones from the nose 
and tail of the configuration do not contribute to the interference wave 
drag. 


Wing and tail-surface components of equivalent bodies of revolution 
have coincident axes (d = 0)j therefore, the dimensionless zero-lift 
interference wave drag between such equivalent -body components can be 
calculated by the present method with values of l(S-y,0). The zero-lift 
interference wave drag between the actual components must, however, be 
obtained from a graphical averaging process. Values of l(&ijj0^ are ~ 
calculated directly from equation (B2) with B = 0. Values of l(ajj,B) 
have been calculated for 0 < a^j < 3.10 in increments of 0.01 and 
0 < B < 3.0 in intervals of 0.2 and are presented in table IV. A value 
of aij =3.1 corresponds to a longitudinal separation of component cen- 
ters of 1.1 half -component lengths of the largest component and should 
be adequate for most practical configurations. 
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APPENDIX C 

AN ANALYTICAL APPROXIMATION OF AIRFOIL-SECTION AND 
PLAN-FORM THICKNESS DISTRIBUTIONS 


The area distribution of a wing or tail surface depends upon the 
thickness distribution of the plan form which , in turn, depends upon the 
thickness distribution of the airfoil section. The plan-form area and 
thickness distributions can not be expressed in analytical form unless 
such an expression is available for the airfoil-section thickness distri- 
bution. In some cases the airfoil-section thickness distribution is 
available in analytical form; for example, that for the NACA four-digit, - 
series airfoil sections is given by a power series. In many other cases, 
such as the NACA 6-series airfoil sections, the thickness distribution is 
not available in analytical form; however, the thickness distribution is, 
in general, accurately specified at many airfoil-section chord stations. 
The slope of the airfoil-section thickness distribution is always known 
for at least one chord station, usually at the chord station for maximum 
thickness, where the slope of the thickness distribution is zero. In this 
appendix a simple approximation of the airfoil-section thickness distri- 
bution will be made which affords an analytical expression for the thick- 
ness distribution from a knowledge of the thickness at specified chord 
stations and the slope of the thickness distribution for at least one 
chord station. With this information, an analytical expression for the 
plan-form thickness distribution of wings and tail surfaces will also be 
determined. The distributions will be determined in the dimensionless 
coordinate system described in appendix A. 


Airfoil -Section Thickness Distribution 


The thickness distribution of an airfoil section is generally 
specified at a sufficient number of chord stations that it can be approx- 
imated quite accurately by a series of adjacent parabolic-arc segments of 
continuous slope. Concerning zero-lift wave-drag calculations, errors 
introduced into the calculation of §(£)> S 1 ( | ) ^ or S"(|) by this approx- 
imation are far less significanc than those inherent in the various 
numerical calculation methods for the wave drag. 
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Consider the dimensionless airfoil-section thickness distribution 
and its derivatives shown in sketch (d). The thickness distribution is 



r'(a) 




Sketch (d) 


expressed analytically as 

t(cl) = f(a) for -1 < a < 1 (Cl) 

where f(a) is a function of a which defines the shape of the thickness 
distribution and a is the dimensionless chord station measured positive 
and negative in the free -stream direction from the 50-percent chord as the 
origin. It should be recalled that the chord station defined in this 
manner is made dimensionless with respect to the half chord of the airfoil 
section. The function f(a) can be determined analytically as a series 
of adjacent parabolic -arc segments of continuous slope when the following 
conditions are satisfied: 

i. Each parabolic segment satisfies the given thickness at two 
consecutive chord stations, a-j_ and a.^ +1 . 



k2 


ii. The slope of the thickness distribution is continuous and is 
known for at least one chord station. 

The function f(a) can be written for a parabolic arc through two 
consecutive points, and a,j_ +1 , such that 

t ( a) = a(a-ai) 2 +b(a-ai +1 ) 2 +c(a-ai)(a-ap +1 ) for a£ < a < aj. +1 

(C2) 

also 

r'(a) = 2a(a-ai)+2b(a-ai+i)+c[2a-(ai+ai+i)] for ai <. a < af+i 

(C3) 



Since -r'(a) is linear and T"(a) is constant over each segment of 
parabolic arc, it is more convenient to express T'(a) at each ai as 

T'(ai) = t' (a i+ 1 )-(a i+ 1 -ai)T"(a) (C8) 

At each predetermined chord station, of the actual airfoil section 

t (a) = T(ai) where r(ai) is known, T'(a) ~ T'(ai) where T'(ai) is given 
by equation (C8), but is known for at least one value of ai, and 
t" ( a) ~ t" ( a^ ) where T M (ai) is given by equation (C4). 

To illustrate the procedure for finding -r(a), r'(a), and T n (a) 
consider a simple parabolic-arc airfoil section and an MCA 6XA00Y airfoil 
section. For the parabolic-arc airfoil section the known conditions are 

i. r(a±) = t(- 1) =0 
ii. T ( a i+1 ) = t(0) = 1 
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iii . t' (a i+1 ) = t* (0) =0 

iv. Oi+i-Cti = 1 

From conditions i and iv in equation (C6) b = 0; from conditions ii and iv 
in equation (C5) a = 1; and from conditions iii and iv and the fact that 
a = 1 in equation (C7) c = -2; hence, from equations (C2) to (C4) 


T(a) = 1-a 2 
r ' (a) = -2a 
r"(a) = -2 


For the NACA 6XA00Y airfoil section the known conditions are 

i. 0.6 < a< 1 T(a) is given for a,i +1 -a^ = 0.1 

T(a) is linear, hence t'(0.6) can he found 

ii. -0.8 ^ a <; 0.6 r(a) given for Oj_ + 1 -a i = 0.1 

iii. -0.95 < a< -0.8 -r(a) given for a,i + 1 -aj[ = 0.05 

iv. -1 < a < -0.95 r(a) given for a^ + 1 -a^ = 0.025 
From condition i, T(aj_) is given and 


t 1 (a*) = - 

T"(ai) = 0 


( a i+i)~ T ( a 'i) 

a i+i“ a i 


(09) 


To continue from condition i, t'(0.6) must he found and used with 
condition ii. From equation (C9) 


t 1 (0.6) = T ( 1 )- T (0- , 6 ) = 2.5[t(1)-t(0.6)] (CIO) 

1 - 0.6 

At a = 0.5, from condition ii and the result of equation (CIO), a, h, 
and c are found from equations (C5) to (07). In equation (C7) 

T > (oi +1 ) = t’ ( 0.6). Next t"(0.5) is found from equation (C4). Then suffi- 
cient information is available to calculate t'(0.5) from equation (C8). 
With t*(0.5) the above steps are repeated and conditions ii to iv are 
used until t(-1), t*(-1), and r"(-l) are determined. 



Wing Plan-Form Thickness Distribution 


kk 


Consider those wings which satisfy the following conditions. 

i. The boundaries of the wing are such that the plan form and the 
spanwise variation of the local chord can be expressed analytically. 

ii. The wing has the same basic airfoil section at all spanwise 
stations; that is, only spanwise variations in thickness and/or thickness 
ratio are allowed. 

iii. The thickness distribution of the basic airfoil section can be 
expressed analytically as described in the preceding part of this appendix. 

iv. The spanwise variation of the thickness along lines of constant 
percent chord can be expressed analytically. 

Actually, condition ii can be relaxed to include wings with one basic 
airfoil section out to a discontinuity in plan form, for example, a fence 
or extended leading edge, and a different basic airfoil section from the 
discontinuity out to the wing tip. In such cases each portion of the plan 
form is treated as a separate entity. 

Consider the one wing panel of a typical plan form shown by the solid 
lines in sketch (e) in dimensionless coordinates. If the abcissa of any 


v 



point in the plan form is measured from the local 50 -percent chord station 
and this distance is divided by the half chord of the local airfoil sec- 
tion, any point in the plan form, (§,tj), can be expressed as a chord 
station of the local dimensionless airfoil section. With the aid of 
sketch (e) this is accomplished by the transformation 


a<U,r)) 


£-KqT) 

5 ( 11 ) 


(Cll) 
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where K 0 is the tangent of the sweep angle of the 50-percent chord line 
of the extended wing panel and d(r|) is the local half chord. All points 
in the plan form which yield the same value of a(£,T]) define lines of 
constant percent chord (see long dashed line in sketch (e.)), and since 
the basic airfoil section is the same at all spanwise stations, the dimen- 
sionless thickness of the local dimensionless airfoil section is constant 
along these lines. Therefore, the thickness at any point in the plan 
form can be found as the product of the thickness at the proper chord 
station of the basic airfoil section and a term which allows for the span- 
wise variation of the thickness along lines of constant percent chord. 

If $(t)) represents the known analytical spanwise variation of the thick- 
ness along lines of constant percent chord, an analytical expression for 
the wing plan-form thickness distribution is 

t(|,T]) = $>(Ti)T[a(!,T])] (C12) 

where T[a(|,T])] is a known analytical expression for the basic airfoil 
section and a(|,T)) is given by equation (Cll). 

Equation (C12) is applicable to all generalized wings which satisfy 
conditions i to iv above. For instance, if the plan form is such that 
the 50 -per cent chord line is nonlinear, equation (C12) is applicable when 
equation (Cll) takes the more general form 

, . S-Kol+fU) 


In this case f(q) describes the nonlinearity of the 50-percent chord 
line. 


Equations (Cll) and (C12) can be modified in a very simple manner 
to yield the plan-form thickness distribution, and hence 

§ (s, e 0) 


and 



for sheared wing components of the equivalent bodies of revolution required 
when calculating the zero-lift wave drag of an airplane -type configuration. 
Consider the complete wing shown in figure 5. The traces in the 

plane of the intersection of the oblique planes and the plane 
define an angle \jr with the rj axis such that 



1+6 




tan 1 f3 


Co 


cos 0 


Hence, the traces through the extended wing tips and the | axis (see 
lines XX and YY in fig. 5) define a longitudinal distance | p(b 0 /c 0 )cos 0 
In accordance with reference 5 each element of the upper (positive 
t] direction) wing panel is sheared a distance |p(b 0 /c 0 )cos ©|tj forward 
and each element of the lower (negative T) direction) wing panel is 
sheared a distance |p(b 0 /c 0 )cos 0 1 r) rearward for 0 <, 0 < n/2 (see the 
dashed lines in fig. 5). From figure 5 it can he seen that the tangent 
of the sweep angle of the sheared 50-percent chord line of the upper and 
lower extended wing panels can he defined, respectively, as 


K u 

= K 0 -3 $2. cos 0 
u o 

(C13) 

K L 

= K 0 +p — COS G 

c o 

(C14) 


The quantity K, rather than 3 and 0, can be used to define the shearing 
of the wing panels . The upper and lower wing panels are related in a 
simple manner hy 


K U+ K L = 2K 0 


(C15) 


It will he recalled that the tangent of the sweep angle of the 50-percent 
chord of the extended wing panel is used in equation (Cll) to define the 
parameter a. Hence, if K 0 in equation (Cll) is replaced hy K, the 
plan-form thickness distribution of the wing panels of the various sheared 
configurations can he written 

t(|,T),K) = ®(ri)T[a(5,Tl,K)] (Cl6) 


where 

a(S,Ti,K) = (017) 

c{r\) 

and K is defined by equations (C13) and (Cl4). 

In calculating the area distribution and derivatives of the area 
distribution for the various sheared configurations, the above means of 
identifying the sheared configurations permits the use of only one vari- 
able, K, to define these distributions. This will be demonstrated by the 
sample calculations of appendix D. 
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The equations derived above for the dimensionless wing plan-form 
thickness distributions are also applicable to tail surfaces when the 
characteristic lengths c 0 /2, b 0 /2, and t 0 /2 for the tail surfaces are 
used to obtain dimensionless coordinates in the x, y , and z directions, 
respectively. For horizontal tail surfaces, equations (C13) and (Cl4) 
are used with = Qy. For vertical tail surfaces, equation (Cl4) is 

used with SyT = 0w+ rt / 2 - The sheared configurations for vertical tail 
surfaces are symmetrical in jt rather than jt/2. 
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APPENDIX D 


AN APPROXIMATE ANALYTICAL METHOD FOR CALCULATING 
AREA DISTRIBUTIONS AND THEIR DERIVATIVES FOR 
SHEARED WINGS OR TAIL SURFACES 


The analytical expressions which will be developed, in this appendix 
for calculating the area distribution, the first derivative of the area 
distribution, and the second derivative of the area distribution are 
applicable to wings or tail surfaces which satisfy the following 
conditions : 

i. The spanwise variation of the 50-percent chord line is known 
analytically. 

ii. The boundaries of the plan form, the spanwise variation of the 
local chord, and that of the thickness along lines of constant percent 
chord can be expressed analytically. 

iii. The same airfoil section is employed at all spanwise stations. 

iv. The thickness of the airfoil section is known at a sufficient 

number of chord stations to permit an analytical approximation of the 
airfoil-section thickness distribution in the manner described in 
appendix C. 

Consider the typical half plan form for wings and tail surfaces 
shown in dimensionless coordinates in figure 6. The figure represents 
all trapezoidal -type half plan forms with taper ratio, 0 < A < 1. The 
general form of the equations which will be developed for the various 
distributions is also applicable to the case for A = 1. The details 
of the solutions for this case, however, are simpler than those for 
0 A < 1. The differences will be pointed out later when some sample 

calculations are presented. The solid lines in figure 6 represent any 

half plan form with 0 < A < 1 , and the short dashed lines represent the 
extended half plan form. It should be noted the actual half plan form 
and the extended half plan form are identical for A = 0. It should also 
be noted that all equations are based upon extended half plan form. 


EXACT EQUATIONS FOR THE DISTRIBUTIONS 


The dimensionless area distribution of any half plan form and its 
sheared configurations is given by 
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r \U,K) 

S(6,K) = / r(£,Ti,K)dTi (Dl) 


Differentiation of equation (Dl) with respect to 5 yields the first 
derivative of the area distribution of the half plan form 


§'(£,K) = 


,T1 u ( ^ K) St(£,tj,K) 


\U>K) 




dT l+ T U> VK) 


- t(|,t 1 ,K) ii 
vs ' *i* at 


(D2) 


Differentiation of equation (D2) with respect to | yields the second 
derivative of area distribution of the half plan form 


i"(S,K) = f 


nJSjK) 




t(i>v k) 


t(6,ti 2 ,K) 


d ^(S<yK) ^ d T ( t^T) u ,K) 

hi 2 11 + ci£ dt dt 

d- g n u _ ^(t,r| z ,K) dr) z d-r(t,T) z ,K) dT] z 
d | 2 ’ St dt dt dt 

d% 

dt 2 


Ik 

dt 


+ 


(D3) 


In equations (Dl) to (D 3 ) t(£,t],K) is the plan-form thickness distribution 
of any of the sheared configurations and is given by 


t(|,T],K) = 0( ^)t[cc( ] 


(C16) 


where 0 (t]) is the known analytical spanwise variation of the thickness 
along lines of constant -percent chord, t[<x(|,t],K)] is the local airfoil- 
section dimensionless thickness at a given percent chord, a(|,rj,K) defines 
lines of constant -percent chord for any point (£,r|) in the plan form and 
is given by 

= f^y (C17) 


and K is given by equation (C13) or (Cl4). The upper and lower limits 
of integration, ^(f^K) and t)^(£,K), respectively, are the boundaries of 

the half plan form. Also 5 t(£;,t],K)/c)£ is continuous between limits of 
integration. 

To evaluate equations (Dl) to (D3) it is convenient to make a change 
of variable from q to a in accordance with equation (C17). In this 
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manner one can subsequently make use of the fact that t"(cx,) is constant 
over each interval of parabolic arc which approximates the airfoil-section 
thickness distribution (see appendix C). With this change of variable 
let 


a = a,(|, t),K) 


(DU) 


\ - V 5 ’ K) 


(D5) 


and from equation (C17) the new lower and upper limits of integration 
become, respectively, 


a >2 = a-U^jK) 
= cx ( 


(d6) 


The change of variable for each term of equations (Dl) to (D3) which 
contains a function of t(|,t),K) can be made by means of equations (Cl6), 
(C17), and (D4) and various combinations of their partial and total 
derivatives. Thus, 


t(S,t),K) 

St(|,t),K) 

hi 

S 2 r(a,ri,K) 

St 2 

dT(6,n,K) 


$(Tj)r(a) 
0(q)T' (a) ~ 


3>('n) T "(a) 

<E>( Tj)x* (a) 



T (o)+0 ( T) )t * ( a) 


eta, dr] 

Sq d| 


(D7) 

(D6) 

(D9) 

(DIO) 


Exact Equation for S(£,K) 


Substituting equations (D7) in equation (Dl) and making the change 
of variable from T] to a 

§(|,K) = J 4>(a,6,K)r(a) ^ da (Dll) 

a Z 

Grouping the integrand of equation (Dll) as t(<x) and 
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® (a, £,K) da 
da 

integration by parts yields 

a n a U 

SU,K) = [r(a)t(a,^K)] a " - J t ' (a)t(a, g,K)da (D12) 

J(X z 

where 

\|f(a,^K) = f »(a,|,K) ^ da (D13) 

J da 

After integration by parts of the second term in equation (D12) 

(Xu o*^u 

g(|,K) = [T(a)i|f(a,5,K)-T'(a)r(a,|,K)] + / T "(a)r(a,|,K)da 

l 'Irv , 

(DlU) 


a Z 


where 


r(a,i,K) = j' i|f(a,5,K)da 
Exact Equation for S'(£,K) 


(D15) 


Substituting equation (D8) and equation (D7), evaluated at q u and , 

in equation (D2) and making the change of variable from T) to a yields 




S'U,K) = J ®(a,6,K)T'(a) || ^ da+d>( H u )t(o u ) djy- - o( n z )r(a z ) 

(Dl6) 

Grouping the integrand of equation (Dl6) as T'(a) and 

$(a ^ jK) If S da 

integration by parts yields 

am 

S'U,K) = [T*(a)fl(a,i,K)]^ - J r"(a)fi(a, ^da+^r^Mau) dk - 


$(q z )r(a z ) 

where 

fl(a,|,K) = f ®(a,6,K) 


d7 l l 

(D17) 

d| 

Sa d 7 , 

— — da 

to 

(D18) 
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Exact Equation for S"(£,K) 


Substituting equations (D9), and equations (D7), (D8), and (DIO), 
evaluated at and in equation (D3), making the change of variable 
from i) to a and collecting the terms 


<x u 

S"(6,K) T"(a)A(a,{,K)a a + ( 2 ^ + ^^i)t(^) T .( CItt )^l + 


■^(^u) far^X d-Tiu 


dT lu V d £ 




' u ' as 2 


(a. ) - ( 2 , *21 d1 h 

(au) t a £ ST IF 


&T)., 


" d< &(11 z ) 

f^l\ 2 




d£ 2 


r(a,) (D19) 


where 


A(a,£,K) 


®(a,|,K) 


Sa\ djrj 

da 


(D20) 


Limits of Integration 


The limits of integration in the equations for the various distribu- 
tions depend upon the plan form being considered. For practical half plan 
forms only four limits of integration need be considered - the chord on 
the £ axis , the leading edge, the trailing edge, and the finite tip for 
^ a ^f plan forms with 0 < A ^ 1. It should be noted that integration over 
the half plan form in the rj direction is transformed, by equation (C17), 
to an integration over the basic airfoil section in the a (or 5 ) direc- 
tion. Furthermore, the limits of a integration never exceed il regard- 
less of the value of £ . The limits of a integration as determined from 
figure 6 and equation (C17) for the actual and extended half plan forms 
are tabulated below: 
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Limits of a integration 

Limit, 

location 

Value 

Type of limit 



Where applicable 

^ axis 

i 

1 

Lower 

-1 < i < 1 

All K 

Leading 

edge 

-1 

Upper 

-1 < | < K 

K > -1 

Lower 

K < i < -1 

K < -1 

Finite tip 
(0 < A < 1) 

!-KT) t 

Upper (actual 
half plan form) 

(Kr^-A) < | < (Kr^+A) 

All K 

s(n) 

Lower 

(extended tip) 

Trailing 

edge 

1 

Upper 

K < i < 1 

K < 1 

Lower 

1 < £ < K 

K > 1 


Presented, in figures 7(a) and 7(t>) for half plan forms with 0 < A < 1 
and A = 0, respectively, are the regions and limits of a integration 
for all possible values of K. It should be noted that the distributions 
of the actual half plan forms with 0 < A < 1 (fig. 7(a)) are found by 
subtracting from the distributions of the extended half plan forms, for 
values of (Kr^-A) < £ < (Kt^+A), the distributions for the extended tip 
defined by tj^ < tj < 1 # It is for this reason that the half plan forms 
are made dimensionless with respect to the extended half plan form. Note 
also that the upper limit of integration is always constant (either -1 
or 1) and the equations for the distributions can be expressed as a func- 
tion of the lower limit. More important, however, this treatment of the 
limits of integration permits the calculation of the distributions for 
the half plan forms with 0 < A < 1 from a very simple manipulation of the 
basic data for the extended half plan forms. This procedure will be 
explained in detail in the sample calculations. Similarly, in the case 
of an airplane-type configuration, the distributions of the exposed half 
plan forms can be conveniently found from the basic data of the extended 
half plan forms . 

It now remains to investigate the effect of the limits of integration 
on the terms in equations for the distributions (eqs. (Dl4), (D17)> and 
(D19)) which are evaluated at the limits. Only the terms in these equa- 
tions which contain t(cx) and/or drj/ dg , evaluated at the limits, will 
vanish when the values of £ are such that the limits of integration of 
the half plan form are 

i. Along the £ axis, in which case di^/dg = 0. 

ii. Along the leading edge, in which case t(o u ) = t(- 1) = 0 or 
t(o, Z ) = T ( -1 ) = 0. 
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iii. Along the finite tip, in which case di^/d^ = 0 or dru/d£ = 0 
when T) t is the lower limit of the extended half plan form. 

iv. Along the trailing edge if the basic airfoil section has or 
is assumed to have zero thickness at the trailing edge, in which case 
t(o u ) = t(1) = 0 or t(cli) = t(1) = 0. 


APPROXIMATE EQUATIONS FOR THE DISTRIBUTIONS 


The approximate equations for S( £,K) , §• ( £,K) , and S" ( £ ,K) will be 
determined by examining equations (Dl4), (D17), and (D19), respectively, 
in the light of the approximation of the airfoil -section thickness dis- 
tribution made in appendix G. It will be noted that each of equa- 
tions (Dl4), (D17), and (D19) contains a term of the form 


a u 

T"(a)f(a,!,K)da (D21) 

l 

For the approximation made in appendix C, t" (cl) ~ t" ( cc.j) = a constant for 
a ± < a < a i+1 (see sketch (d)). Therefore, if the integration in expres- 
sion (D21) is performed as the sum of integrations over the intervals 
a i+i“ a i> T " ( a i ) can be taken outside the integral sign for each interval. 
The sum of the individual integrations which make up the integration of 
the expression (D2l) can have two forms, depending on whether a u is 
less than or greater than <x z . The form of the integration for these 
two cases can be determined with the aid of figure 8 which defines the 
stepwise limits and direction of a integration. Case I, for which 
o-u < a z applies, regardless of plan form, when £ < K; and Case II, for 
which a u > a ^ applies, regardless of plan form, when £ > K. 

With the aid of figure 8 and the approximation in appendix C, the 
expression (D2l) can be written 



(Xu 

J r"(a)f(a,|,K)da * T "(a z ) f 

CLn nr.-, 


a i(a z ) 


U+i 


i=i(a z ) 


f(a,£,K)da + 
a i-i 

"(a-i-i) J f(a,|,K)da + 


°± 


C a u 

T n (au_x) / f(a,^,K)da for Case I 

Jq u 


(D22) 


and 
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CX/Q i 

J r"(a)f (a,£,K)da x t "(< x 2 ) J 
a 7 a 7 


a i(a^ )+i 


U-i 


f(a,S,K)dot + 
, a i+i 


^ T"(ai ) J f(a,S,K)da 
i=i (aj )+i 


r"(a,y) / f(a, £,K)da for Case II 
a U 


(D23) 


If aj and can be made to coincide with and ay for Case I and 

with <Xj_ +1 and ay for Case II, the first and last terms of equations (D22) 
and (D23) will vanish. This can "be accomplished for the first terns by 
choosing £ such that cc^ always corresponds to a value of aj_ for 
which the airfoil-section thickness is specified, and for the last terms 
by always integrating to the leading or trailing of the extended half 
plan form (i.e., = +l). As mentioned earlier, this is accomplished 

for half plan forms with 0 < X < 1 by subtracting the distributions of 
the extended tip from those for the extended half plan forms. 


The approximate equations for S(g,K), S'(£,K), and §"((-, K) are 
found by substituting the form of equations (D22) and (D 23 ), with the 
proper integrand, in the integral term of the exact equation for each 
distribution. For §(£,K), equation (Dl4), f(a,£,K) = r(a,g,K). For 
S'(g,K), equation (D17), f(a,5,K) = ft(a,S,K). For 3"(£,K), equation (D19), 
f(a,£,K) = A(a,6,K). The terms of equations (D14), (D17), and (D19) which 
are evaluated at the limits of integration are treated in the same manner 
as discussed in connection with the exact equations for the distributions. 


SAMPLE CALCULATIONS 
Complete Half Plan Forms, 0 <* X < 1 


The area distribution, the first derivative of the area distribution, 
and the second derivative of the area distribution will be calculated for 
the families of half plan forms shown in figure 7. The airfoil sections 
will be considered to have a linear spanwise variation in thickness along 
lines of constant percent chord, and, for simplicity, will be considered 
to have zero thickness at the trailing edge. 3 The spanwise variation of 
the half chord is given by 

generally, the thickness at the trailing edge of the airfoil section 
is sufficiently small that it can be considered to be zero without intro- 
ducing any appreciable error in the distributions. Any small error which 
might be introduced, however, has a negligible effect upon the zero-lift 
wave -drag calculations. 



56 


c(n) = i-n 

(D24) 

From equations (C17) and (D24) 


a . ^ 

(D25) 

1-T) 


and 


i £- K 

1 = 1 - ^ic 

(D26) 


The linear spanwise variation of the thickness along lines of constant 
percent chord is 


®(tj) = 1-CT) for 0 < r\ < 1 (D27) 

where 

C = i=2. for 0 < A < 1 (D28) 

1-A 

In equation (D28) cr is the ratio of the airfoil-section thickness at the 
finite tip to that at the vertical plane of symmetry, and A is the taper 
ratio of the actual half plan form "being considered. For the spanwise 
thickness variation, C is normally defined by only the numerator of equa- 
tion (B28). In ^e present case, however, the factor l/l-A is required 
in order that equation (D27) he valid for 0 < rj < 1. This factor is the 
ratio of the extended semispan, b 0 /2, to that of the actual semispan, h/2. 
With C defined as in equation (D28), for half plan forms with 0 < A < 1, 
all or a portion of the spanwise thickness distribution of the extended 
half plan form may be negative. This imposes no serious consequences 
since the distributions for 0 < A < 1 are obtained by subtracting the 
contributions for rj^ < r\ < 1 from the distributions for 0 < r\ < 1. 

For extended half plan forms it is always possible to calculate the 
distribution at values of £ for which the lower limit of integration 
corresponds to values of ccj_ for which the airfoil-section thickness is 
specified, and the upper limit of integration is always the leading or 
trailing edge. Hence, for extended half plan forms the first and last 
terms of equations (D22) and (D23) vanish. With this information plus 
conditions i to iv previously discussed in connection with the limits of 
integration, the approximate equations from which the distributions can 
be calculated for half plan forms with 0 < A < 1 are 
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a u 


S(|,K) « [T(a)\|r(a, £,K)-t' (a)r(a, £,K) ] + 


r 2 


a i-i 


T"(ai_!) I r(a, |,K)da for Case I 
i=i(a z ) ai 

1 „ a i+i 

r"(ai) / r(a, |,K)da for Case II 


> 0)29 ) 


V,i=i(a z ) +1 


a n * 


S f (^K) » [T f (a)fl(a,g,K)] 


a u 

a. 


r 2 


ai _ 


^ T "( a i-i) J H(a ; ^K)da 

i=i(a z ) a i 




for Case I 


> (D30) 


I-i 


for Case I 


§"U,K) «< 


2^ T"(cti) J fi(a,£,K)da for Case II 

v s i=i(a z )+i ai 

C 2 a l-i 

^ T , '(ai_ 1 ) J A(a,5,K)da 

i=i(a 2 ) a i 

I-i a i+i 

/ T"(ai) / A(a,£,K)da for Case II 

^i=i(a z ) + x 


dt|, d| 


d£ 


(D31) 


Sample calculation of S(g > K) .~ For half plan forms with 0 < A < 1, 
first consider the extended half plan form and make a change of variable 
in equation (D27) by using equation (D 26 ), hence 


®( a,f,K) = (1-C)+C 

CC-iV 


<D32) 


From a differentiation of equation ^D26) with 1 respect" to a. 
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Substituting equations (D32) 
indicated integration 

|,K) 


in = i- K 

da “ (a-K) 2 

and (D33) in (D13) and performing the 



(D33) 


(D3^) 


Substituting 

integration 


equation (D34) in equation (D15) and performing the indicated 


r(a,6,K) = ~(l _ C ) ( | -K) In | a-K | + | 


(D35) 


Substituting equation (D35) in the second term (for Case I) of 
equation (D29) 


^ T ,, (a i _ 1 ) J h(a, | ,K)da 

i=i(a 1 ) 




2 

= -(l-C)(|-K) ^ t" (aq ( In |ap_ 1 -K | -1 ) ■ 


i=i(a z ) 

2 

(ai-K)(ln|ai-K|-l)] + | (|-K) 2 ^ T n (ai-1) (Zn|ai_!-K| -Zn|ai-K| ) 

i=i(a i ) 

for Case I (D36) 


Each of the summations in equation (D36) is very large for a close to 
the value of K and infinite for a z = K. To avoid this singularity, 
it is therefore desirable to multiply and divide both terms of the right- 
hand side of equation (D36) by a,j-K. Furthermore, since the uppermost 
limit of the stepwise integration (ay) is -1 along the leading edge of 
the half plan form, equation (D36) can be written as a function of only 
the lower limit, and K. Thus 

)Bj(o, 2,K) for Case I 


2 

I 


,01-1 


t" (cti_i) 


i=i(po 7 ) 


OM 


f(a, | ,K)da = -(l-C) ( 


l-K 

2 W-K 


(D37) 
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where 


and 


A I (a l ,K) = (a z -K) £ t" [ (a^-K) ( Zn l^^-K | -l) - 

i=i(a z ) 

(ai-K)(Zn |ai-K| -1)] for Case I 


Bj(a z ,K) = (a z -K ) 2 ^ t" (a^J ( Zn | -Zn |ai-K | ) 


(D38) 


i=i K) for Case I (D39) 

Similarly, since a y = 1 along the trailing edge, 

I — i a^ +1 . 

£ T ,, (a 1 ) J r(a-,|,K)da = -(1-C) ^jAnfa^K) + 

i=i(a i )+i 0,1 . v 2 

§ f“V) BlI ( a Z; K ) for Case 11 

(D4o) 


where 


I-i 


Ajj(<x^jK) = ( cx, ^ -K ) r"(a 1 ) [ | l) “ 

i=i(a^ ) + i 

(aj[-K)(Zn |ai-K| -1) ] for Case II 


(d4i) 


and 

I-i 

Bn(a z , K ) = (a z -K f ^ t" (c^ )( Zn |a i+1 -K | -Znlc^-K | ) 

i=i(a^)+ 1 f or Case II (D42) 

Substituting equations (D34), (D35), (D37 )j and (D40) in equation (D29) 
an expression for the area distribution of the extended half plan form 
is obtained as 

S(g,K) « (1-C) [C(a z ,K)-D(a z ,K)] + | (“^) t E ( a Z ^ K ) +F ( a Z > K ) ] 

(1*3) 



where 


C(aj,K) = r(a 2 ) - T(a u )-A(a z ,K) (D44) 

D (aq>K) = (a z -K)[r' (a 2 )Zn|a z -K| -t* (a 11 )Zn|a tl -K| ] (D45) 

E(a z ,K) = r(a z ) - t(<x u )+B(o. z ,K) (D46) 

F(a z ,K) = (aj-K) T*(a z ) “(“^ T '(“u) OW 


The left -hand side of equations (D44) to (D47) is correctly given as a 
function of only a z and K even though a u , t(o u ), and t'(o/ U ) appear in 
the right-hand side of the equations. For the extended half plan form it 
will be recalled that ct^ = ay is always +1, t^) (included for the 
general case) is zero along the leading edge but may not be zero along 
the trailing edge of the half plan form, and T^dy) is one constant value 
along the leading edge and another constant value along the trailing edge 
of the half plan form. For these sample calculations, it will be recalled 
that t(<x u ) has been assumed zero along the trailing edge; consequently, 
the term containing t(o u ) in equations (D44) and (D46) vanishes. For 
Case I and Case II A(a z ,K) and B(a z ,K) are given by equations (D38) and 
(D39) and equations (D4l) and (D42), respectively. 

To evaluate equation (D43) it will be found convenient to calculate 
the values of A(a z ,K) through F(a z ,K) in tabular form for all values of 
K of interest and for all values of <x z corresponding to values of 
-1 <1 ap 1 for which the basic airfoil-section thickness distribution 
is given. Then, depending upon the value of a z is obtained from the 

table on page 53 and equation (D43) is evaluated with the aid of the 
tabulated values of A(a z ,K) through F(a z ,K). 

Equation (D43) is also applicable to half plan forms with A = 0, in 
which case the actual and extended half plan forms are identical. 

Equation (D43) is applicable to all geometrically similar extended 
half plan forms . One such geometrically similar half plan form is the 
extended tip of the extended half plan form defined by tj+ £ tj <* 1. Hence, 
the area distribution of the extended tip can be easily found if, .In 
equation (D43), a z of . the extended-tip half plan form is taken -o^. 
of the actual half plan form. Then S(i,K) for the extended tip^ecomes 
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f(l-C) 


3(6, K) = { 


[ C (^ > K)-D(a t ,K)] + §(^|) [E(at,K) + F(at,K)] 

for (Kip^-A) < | < (Ktj^+A) > 
Equation (D43) for | ^ (K^±A) 

4 

(D48) 


If happens to correspond to values of ap for which the airfoil- 

section thickness is given, C(clj,K) through F(a z ,K) have previously been 
calculated in finding S(£,K) for the extended half plan form. Otherwise, 
these quantities can be obtained from individual or combined plots of the 
tabulated values of C(a z ,K) through F(a z ,K) versus their arguments; that 
is, the values of these quantities are read from the plots at values of 
0,7 obtained by letting a z = op. Hence, the area distribution of the 
half plan forms with 0 < A < 1 is given by the difference between equa- 
tion (D43) for all values of £ and equation (d 48) for 
(Kt^-A) < 6 ^ (Kt^+A). 

Sample calculations of 5 1 ( £ ,K ) . - Again the extended half plan forms 
will be considered first. All the terms of equation (D18) except 6 a/o£ 
have been found in terms of cc, £, and K in the sample calculation of 
S(|,K). This term is evaluated from equations (D25) and (D26) as 


6 a, _ a-K ( d49) 

8 £ 6 -K 


Substituting equations (D32), (D33)j and (D49) in equation (Dl 8 ) and 
performing the indicated integration yields 


ft(a,|,K) 



(D50) 


From a comparison of equations (D35) and (D50) it can be easily determined 
from an inspection of equations (D 36 ) to (D39) that by substituting equa- 
tion (D50) in the second term, for Case I, of equation (D30) one can obtain 


d 

i=i (cl 1 ) 


CCn 


X -1 


fi(a, £,K)da 


“i 


(1-C) An;(a z ,K)-C 

a z -K- (a z -K) 

for Case I (D5l) 


and similarly 
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1-1 


, a i+i 


^ T'^ai) J fl(a,|,K)da 

i=i(a z )+i 


a i 


(1 " C) ^=K A lf a Z' K ) 
C 7 ~ K 'g B Il( a Z,K) 

(a-Z-K) 2 


for Case II 
(P52) 


where A(aj,K) and B(<xj,K) are given for Cases I and II by equations (D38) 
and (D39) and equations (D4l) and ( D42 ) , respectively. Substituting 
equations (D50), ( D51 ) , and (D52) in equation (D30) an expression for the 
first derivative of the area distribution of the extended half plan form 
is obtained as 


S* (g,K) - -(i-c) 


a Z" K 


[A(a z ,K)+D(az,K)]+C 


l-K 


( a z - k) 


g [B(aj,K)+F(a 2 ,K)] 

(D53) 


As in the case of the area distribution, S f (£,K) for the extended half 
plan form (actual half plan form for A = 0) Is obtained from equa- 
tion (D53) by using the proper value of a 1 obtained from the table on 
page 53 for each value of £ for which § ! (£,K) is calculated. Values 
of A(a^,K), B(a^,K), D(a^,K), and F(clj,K) are obtained from the previous 
calculations of 3(£,K). For the extended-tip half plan form § T (^K) Is 
obtained by following the same procedure outlined in finding S(£,K). 

That is, in the region for which (Kn t -A) < | < (KT] t +A), a l in equa- 
tion (D53) is replaced by a^, the values of A(a^,K), D(a^,K), and 
F(a z ,K) are obtained from the aforementioned tabulated values, or plots, 
by letting and the resulting values of S T (£,K) are subtracted 

from those of the extended half plan form at the proper values of £ . 
Obviously, for £ ^ (Kr^_±A), S*(f-,K) for the extended tip Is given by 
equation ( D53 ) . 


Sample calculation of 3 n (£,K) .~ For the extended half plan forms, 
substituting equations (D 32 ), (D33)> and (D49) in equation (D20) yields 
the function 


A(a,|,K) 


l-c _c_ 

|-K a-K 


(D54) 


Using equation (D54), the first term of equation (D31) becomes 

[t’ (a u )-T* (a z )] + B( a 2 >K) (D55) 

where B(aj,K) for Cases I and II is given, respectively, by 
equations (D39) and (D42). 
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To evaluate the second and third terms of equation (D31) consider 
the equations of the leading and trailing edges of the extended half pi an 
forms . In terms of both tj and a variables 

= Mr (D56) 

and 

a = +1 (D57) 

The upper and lower signs of the ± sign in equations (D56) and (D57) are 
used with the leading and trailing edges, respectively. First consider 
the evaluation of the second term of equation (D31) when the leading or 
trailing edge is the upper limit of integration. From equation (d 49), 
evaluated at a u , and equation (D57) 


_ K±1 

W ~ kT 


(D58) 


Since 60^/^^ = da^/dT^ = l^dr^/da^), from equation (D33), evaluated 
at a u , and equation (D57) 


6a u _ (K±l) 2 
6t) u " K-l 

From a differentiation of equation (D5 6) with respect to | 


(D59) 


d % _ 1 

d| K±1 

From equation (D32), evaluated at a u , and equation (D57) 


(D60) 




(1-C)+C 


K-| 

K±1 


(d 6 i) 


Using equations (D58) to (D6l) one obtains for the second term (upper 
limit) of equation (D3l) 


3t* ( a u ) 

i-K 


( C|-K)+ ( I-C ) 
K±1 


(D62) 


A general expression identical to that of (D62) is similarly obtained for 
the last tern (lower limit) of equation (D3l). 


For the extended half plan forms S"(|,K) is found by properly 
combining expressions (D55) and (D62) (and its counterpart for the lower 
limit) for certain ranges of the value of |. Regardless of the value 
of K, for -1 ^ | 1 , the last term of equation (D3l) vanishes, because 

dT] z /d| = 0 for these values of |. For | K | > 1, and for ||| >1, when 
the leading edge of the extended half plan form is the lower limit, the 
trailing edge is the upper limit of integration and vice versa. There- 
fore, for the extended half plan form (actual half plan form A = 0) 
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S"U,K) ={ 


~ (Kt2)(1-C)+3(CS-K)' 

(K±l)(6-K) 


'(KT2)(1-C)+3(C&-K)' 


T ' (a ^ )+ ^K B(a ^ K) ‘ T ' (a ^ 

for -1 £ | £ 1 


(D63a) 


(K±l)(|-K) 

(K±2) (l-C )+3 (C| ~k) 
(KTl)(|-K) 




T ' ( a 'u)+ B(a l 7 K) 


''( a z) 


for M > 1 


(D63b) 


where the upper of the double signs is employed when the leading edge of 
the extended half plan form is the upper limit of integration and the 
lower sign is used when the trailing edge is the upper limit. It should 
be noted that equation (D63b) is constant for C = 1 (constant airfoil- 
section thickness ratio in the spanwise direction). 

To calculate S"(|,K) for the extended tip half plan form, the same 
procedure described in finding §(£,K) and S’(£,K) is employed. Again, 
a i is replaced by for the actual half plan form in equation (D 63 a) 

for values of (Kt^-A) < £ <1 (Kr^+A) and for § ^ (Kr^tA) equation (D 63 b) 

is used. The resulting S"(£,K) is subtracted from that of the extended 
half plan form. 

If the half plan forms considered previously have a more general 
airfoil section such that the airfoil-section thickness is not zero at 
the trailing edge, equations (D44) and (D46) are valid as shown, but for 
the values of £ for which the trailing edge of the half plan forms is 

a limit of integration, the terms of equations (D 17 ) and (D19) which were 

omitted in the above calculations must be evaluated. 



K 

1 

/ , 

r, * fl 

/ 


Tronsformotion = 

/ 

a = f-K77 

/ 



span, b/2. Hence 
as shown in sketch 


Complete Half Plan Forms , A = 1 


For the special case of half 
plan forms with A = 1, the extended 
half plan form is semi -infinite and 
a different definition of K and 
different equations for the bound- 
aries of the half plan form are 
needed. As noted in appendix A, the 
spanwise coordinate of half plan 
Sketch (f) forms with A = 1 is made dimension- 

less with respect to the actual semi- 
the dimensionless half plan form for A = 1 is defined 
(f). With the half plan form so defined. 
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C = 1-cj for A = 1 

and the equations for the distributions derived above are, in general, 
applicable for A = 1. However, for values of £ in the region of the 
finite tip, (K-l) <> £ (K+l), the value of the upper limit defined by 

the a transformation noted in sketch (f) must be used in equations (D44) 
to (D47), that is, 


= 6-Kt^. = 6-K 


Furthermore, unless the values of £ 
at corresponds to the values of a-j_ 
thickness is given, the last term of 
vanish and must be evaluated. 


in this region are so chosen that 
for which the airfoil-section 
equations (D22) and (D23) will not 


Exposed Half Plan Forms 


In calculating the distributions of the exposed half plan forms of 
airplane -type configurations it is convenient to first calculate the dis- 
tributions of the extended half plan forms. This provides an accurate 
determination of plots of A(a^,K) through F(a^ ,k). if % = % d) 
defines the body plan form, integrals of the type 


^ u ~ +1 

CXq =a( T)-g ) 


for the exposed portion of the extended half plan form can be evaluated 
in the manner described above for calculating the distributions of the 
half plan forms with 0 < A < 1. For values of £ at which the distribu- 
tions may be desired, it is unlikely that a(T>g) will correspond to values 
of ai for which the airfoil-section thickness is given. Hence, values 
of A(a^,K) through F(a^,K) are obtained from plots of these quantities 
by letting = a(rjg). Furthermore, for values of £ in the region of 
the juncture of the half plan form and body, the terms containing 
drj^/dl which were omitted above in finding S f (£,K) and S n (£,K) (see 
eqs. (D17) and (D19), respectively) must be evaluated since, generally, 
for practical configurations, dr)^/d£ ^ 0. If the local values of dr^/d^ 
are small enough to permit the approximation dr^/d^ = 0, it will be 
found convenient to treat the exposed extended half plan form as a complete 
extended half plan form and use equations (D29) to (D31). In so doing, 
however, care must be exercised in evaluating k x , ky, k z , Sr, and C and 
in properly shifting the upper (positive tj direction) and lower (negative 
T) direction) sheared half plan forms when combining their distributions 
to calculate the zero -lift interference wave drag between the sheared 
half plan forms. 
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The approximate distribution equations in the sample calculations 
were derived by integrating the basic equations with respect to a in 
order to take advantage of the approximation made in appendix C for air- 
foil sections with nonanalytic thickness distributions. Although these 
approximate equations may be used for half plan forms for which the 
airfoil-section thickness distribution is given by a simple analytical 
expression, it is found more convenient, in such cases, to perform the 
integration of the basic distribution equations in terms of rj as indi- 
cated by equations (Dl) to (D3). For example, the thickness distribution 
of a parabolic-arc airfoil section, as shown in appendix C, is 


t(cx) = 1-a 2 for -1 < a < 1 


From equation (C17) 


T[a(t,T},K)] = 1 


Therefore, from equation (D7) 


t(S,tj,K) 


®(ti) 



Hence, with the spanwise variation of the thickness, o(rj), and the 
boundary of the half plan form specified (with the boundaries known, the 
spanwise variation of the half chord, c(t]), is automatically known) all 
the terms in equations (Dl) to (D3) can be evaluated exactly. Furthermore, 
if dri/d£ for the boundary is continuous, as in the case of an elliptic 
plan-form wing, the calculation of the distributions is greatly simplified 
by virtue of the fact that a single equation specifies the limits of 
integration. 
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APPENDIX E 

A PROCEDURE FOR CALCULATING LIFT DISTRIBUTIONS 
FOR SHEARED WINGS 


The lifting -pressure distribution supported by a wing can be defined 
as the sum of several superposed distributions. These superposed lifting- 
pressure distributions can be written in coefficient form as 

N M 

C p (x,y) = ^a n x n + 

n=o m=o 

where x and y are the conventional rectangular coordinates . The constant 
coefficients a n and b m are determined by specified boundary conditions , 
such as lift, pitching moment, and desired properties of the lift distri- 
bution. In reference 9 a method was developed, independent of thickness, 
for warping delta and arrow plan— form wings with subsonic leading edges 
and supersonic trailing edges to support a desired lift distribution at 
0° angle of attack. The method employed a particular form of equation (El) 
which was considered adequate to determine the lifting— pressure distribu- 
tion of this specialized group of wings with acceptable accuracy. In terms 
of the dimensionless coordinates defined in appendix A of this report, 
this particular form of equation (El) can be written 

C p (S,T)) = a Q + g+bi hl+baTi 2 (E2) 

The longitudinal lift distribution is then easily obtained from 

E(S) = / Cp(|, 7 i)dTi (E3) 

Vs) 

where C p (|,T]) is given by equation (E2), and t^U) and r^U) are 
boundaries of the dimensionless plan form. 

Just as in the case of the 
area distribution, the lift dis- 
tribution can be conveniently 
determined in terms of that for 
the sheared wings of each equiv- 
alent body of revolution. The 
local lifting pressure at any 
point in the plan form is inde- 
pendent of the shearing process. 

Therefore, with the aid of 
sketch (g) it is clear that the 


v 
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local pressure at any point in the sheared plan form, K, can be related 
to the same local pressure at the corresponding point in the original 
plan form, K 0 , by the following transformation 

i op = £-(K-K 0 )tj (E4) 

where the subscript op denotes the original plan form. Using equa- 
tion (E4) in equation (E2) the distribution of the lifting -pressure 
coefficient for sheared wings can be written 


C p (|,T),K) = a D + I + 


K-K 0 

K 0 +l 



|ti|+b 2 T i 2 


(E5) 


The integral of equation (E5) over the plan form in the t] direction 
results in the frontal projection of the force distribution intercepted 
on the given wing by the parallel oblique planes. The direction of the 
frontal projection of the net force is perpendicular to the oblique 
planes (see sketch (a)). With the subscript p used to denote this fact 
integration of equation (E5) yields 


LpU,K) 




+ 



-Mt,K) 

Ml, K) 


(e6) 


The values of the limits of integration for the sheared wings are deter- 
mined in the same manner described in appendix D for area distributions 
of sheared wings. However, in equation (e 6) the limits are expressed in 
terms of the tj integration rather than the a integration. Finally, 
it is the component of the force distribution in the £ direction (i.e., 
the lift distribution) which is employed in equation (2) in this case. 
Hence, from sketch (a) it is clear that the lift in the £ direction is 

LU,K) = LpU,K)sin 0 (E7) 

and 

L'(|,K) = Lp'(|,K)sin 0 (E8) 

Upon substitution of equation (E8) in the dimensionless form of 
equation (2) of this report, it is clear that the total wave drag is most 
conveniently computed by separating the effects on the wave drag of 
thickness and lift, because only the lifting effect is affected by sin 0 
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APPENDIX F 


SPECIAL COMPUTING PROCEDURE FOR THE METHOD 
OF REFERENCE 3 


If equation ( 3 ) of this report is integrated by parts, the zero -lift 
wave drag of any slender configuration can be expressed as a function of 
the first derivative of the area distribution of the configuration, 
S'(x,p,0). From an expansion of S'(x,p,0) in a Fourier sine series, 
the infinite series form of the supersonic area rule was developed in 
reference 2 for calculating zero-lift wave drag for slender-body thin- 
wing combinations; that is, 


where 


. 2 it 


D(P) = K 


00 

^ n[A n (p,0) ] 2 d0 


n=i 


A n (p,e) 



6s(x,p,e) 

Sx 


sin(ncp)dq 


(FI) 


(F2) 


and 


x 


z 1 (p,e)+z(p,e) 

2 


cos cp 


(F3) 


The zero -lift wave drag of the equivalent bodies of revolution is given , 
of course, by 

D( Pj 0) -i q yn[A n (p,0 )] 2 (F4) 

n=i 

The use of the Fourier series requires that application of the supersonic 
area rule be restricted to configurations for which each equivalent body 
of revolution satisfies the condition given by equation (4), 

S'(-l x ,£,9) = S* (Z ,3,0) = 0 

In reference 3 the Tchebichef form of harmonic analysis was employed to 
evaluate equation (F2). This permits zero-lift wave-drag calculations 
by the supersonic area rule from a knowledge of the area distribution of 
a configuration, S(x,p,0). 
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The electronic-machine computing equipment at Ames is programed to 
evaluate the Fourier sine series coefficients and compute the total zero- 
lift wave drag of an equivalent body from the values of the area distri- 
bution determined by a telereader at 201 points equally spaced over the 
total length of an equivalent body. Equation (f 4) requires an infinite 
number of terms of the series. However, the practical application of 
equation (F4) requires that the series be terminated at some value of n, 
say N, for which the series will converge. The existing machine- 
computing program is set up to provide zero-lift wave -drag solutions for 
N up to and including 49. In reference 3 it was demonstrated that for 
configurations with no discontinuities in the area distribution or the 
first derivative of the area distribution of a configuration N = 25 
provides adequate convergence of the series, and the method yields accu- 
rate values of total zero-lift wave drag. This was found to be the case 
for pairs of bodies of revolution of this report which also had no singu- 
larities in the second derivative of the area distribution at the ends of 
the bodies. For pairs of bodies, at least one of which had singularities 
in the second derivative of the area distribution at the ends, all of the 
49 terms available were required to obtain adequate convergence of the 
series. 4 For this reason 49 terms of the Fourier series for a 201 -point 
analysis of the area distribution of a configuration were used in this 
report in order to consistently obtain the most accurate results by the 
method of reference 3. 

In addition to the choice of N, there exists another factor which 
affects the convergence of a Fourier sine series. In reference 10 it is 
pointed out that in performing a Fourier series analysis of the derivative 
of the area distribution of a configuration the higher harmonics of the 
series will be suppressed if the smallest allowed value of the length of 
the configuration is used. Suppression of the higher harmonics results 
in a more rapid convergence and a more accurate result for a given N 
of the series. In the case of multiple -component configurations this is 
easily accomplished (as suggested in ref. 10) by analyzing the area dis- 
tribution of each component over its individual length to obtain the 
zero -lift wave drag of each component alone and by analyzing the combined 
area distribution of any two components over their total combined length 
to obtain the zero-lift interference wave drag between components. The 
effectiveness of this type of individual-length analysis has been 
adequately demonstrated in the EVALUATION OF THE METHOD section of this 
report . 

The machine -computing program for the method of reference 3 does not 
yield the zero-lift interference wave drag directly. However, in addition 
to the total zero-lift wave drag, the program output sheets list the 

4 For bodies lxl and 5x5 in the arrangement of figure 3(a), the series 
actually converges to the analytical value in the first two terms. How- 
ever, to be consistent, data for these bodies computed by the method of 
reference 3 are presented for N = 49. 
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Fourier coefficients. The zero-lift interference wave drag can be 
obtained by proper application of these coefficients. For instance, the 
Fourier coefficients for a two-component configuration can be written 


A T n (P> 0 ) = A E n (P> 0 ) +A H n (P> 0 ) 


(F5) 


From substitution of equation (F5) in equation (F4) the zero-lift inter- 
ference wave drag between components of an equivalent body of revolution 


is given by 




00 

= | q^nA En ((B,0)A Hn (|B,0) 


n=i 


(F 6) 


Hence, by means of the Fourier coefficients of the program output sheet, 
the zero-lift interference wave drag can easily be calculated by manual 
computations in accordance with equation (f6). 


The procedure described above was employed to calculate the 
dimensionless zero-lift interference wave drag at M = 1 for pairs of 
bodies of this report with 6=0 (figs. 3( a ) an l 3( c ))- H was found 
convenient to compute these values according to the relation 


&EH = (DEH) n even +( & EH) n odd 


Body E 


(F7) 


Body H 




The separate terms of equation (F7) 
were also used to compute the inter- 
ference at M = 1 for the pairs of 
bodies of this report with 5=2 
(fig. 3(h)) from the relation 

^EH ” (^EH) n even ^^-^n odd 

(F8) 

The derivation of equation (F8) 
can best he understood by consider- 
ing the two bodies shown by the 
solid lines in the upper part of 
sketch (h). If the area distribu- 
tions of bodies E and H are ana- 
lyzed in a Fourier series over the 
length -1 £ | £ 1, the contributions 
to the total area distribution of 
the even and odd parts of the 
Fourier series, S e (£) and S 0 (£), 
respectively, can be constructed for 
each body as shown in the sketch. 


S(£) 



Sketch (h) 


i 
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The zero-lift interference wave drag between bodies E and H consists 
of the sum of the zero-lift interference wave-drag contributions of the 
even parts of the area distribution and that between the odd parts of the 
area distribution (see eq. (F7)). The zero-lift interference wave drag 
between the even and odd parts of the area distribution is identically 
zero (see, e.g., bodies 1x5 of table II ). If body H were positioned 
as shown by the dashed lines (5 = 2), the contribution to the total area 
distribution of body H of the even part of its area distribution would 
be the same as before, but the contribution of the odd part of its area 
distribution would be given by the dashed curve. This dashed curve is 
exactly the negative of §jj (£) for the original position of body H. 
Consequently, the zero-lift interference wave drag at M = 1 for the 
configurations with 8=2 was calculated by using equation (F8) and 
the data previously obtained from equation (F7) for the configurations 
with 8=0 but analyzed over a length equal to the total combined length 
of the bodies with 8=2. Although this procedure has been employed 
herein to obtain M = 1 interference drag for only 8 = 2, it should be 
noted that the general procedure is applicable for all 8 > 0. 
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TABLE I.- ZERO-LIFT WAVE-DRAG EQUATIONS FOR SEVERAL CONFIGURATIONS 


7k 



>r this configuration only, S 1 ( 2 ) f 0. 

3, W^P^ W^P 2f and P^B are treated in the same manner as W X B, WjP^, and P X B, 

respectively, because of the symmetry of the sheared configurations for 0 £ G £ * and 
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TABLE II . - DIMENSIONLESS ZERO-LIFT INTERFERENCE WAVE DRAG AT 
MACH NUMBER 1 FOR CONFIGURATIONS OF FIGURE 3 


1 

2 

3 

4 

5 

6 

7 

8 

9 


or 

<=> — 

<=> 


— — 

— 



Bodies 

5 = o 

5 = 2 

5=0 

5=0 

5 = 2 

5 = 0 

5=0 



Analyl 

bical 

Method of reference 3 

Pr< 

2 sent met hoc 

1 

l*xl* 

1*X2 

l*x3 

l*x4 

l*x5* 

l*x6 

l*x7 

l*x8 

7.064 

5.892 

5.192 

4.636 

0 

0 

0 

0 

1.068 

7.0o5 

6.909 

5.885 

5.153 

4.636 

-.208 

0 

.004 

0 

1.053 

.538 

.359 

.280 

.650 

.245 

.120 

.074 

7.1 
5.9 

5.2 

4.6 

7.1 

5.9 

5.2 
4.6 

0.98 

.51 

.35 

.28 

.59 

.22 

.12 

.074 

2x2 

2x3 

2x4 

2x5* 

2x6 

2x7 

2x8 

7.368 

7.733 

7.744 

0 

0 

0 

0 

.316 

7.368 

7.367 

7-728 

7.742 

.007 

.004 

0 

0 

.318 

.234 

.191 

.286 

.119 

.068 

.045 

7.367 

7.730 

7.741 

7.358 

7-728 

7.736 

.316 

.234 

.191 

.26 

.119 

.069 

.045 

3x3 

3x4 

3x5* 

3x6 

3x7 

3x8 

9.014 

9.658 

0 

0 

0 

0 

.178 

9.014 

9.015 

9.658 

.001 

0 

0 

0 

.178 

.148 

.166 

.081 

.049 

.033 

9.012 

9.650 

9.003 

9.650 

.178 

.148 

.16 

.081 

.048 

.033 

4x4 

4x5* 

4x6 

4x7 

4x8 

10.86 

0 

0 

0 

0 

.123 

10.86 

10.85 

0 

0 

0 

0 

.124 

.123 

.062 

.038 

.026 

10.85 

10.84 

.124 

.12 

.062 

.038 

.026 

5*x5* 

5*x6 

5*x7 

5*x8 

4.710 

2.946 

2.061 

1.545 

-.487 

4.733 

4.549 

2.942 

2.063 

1.546 

-.452 

-.158 

-.070 

-.039 

4.7 

2.9 

2.1 

1.5 

4.7 

2.9 

2.0 

1.5 

-.42 

-.14 

-.067 

-.041 

6x6 

6x7 

6x8 

2.947 

2.578 

2.213 

-.055 

2.947 

2.945 

2.578 

2.213 

-.056 

-.029 

-.018 

2.946 

2.578 

2.212 

2.941 

2.573 

2.209 

-.055 

-.030 

-.018 

7x7 

7x8 

2.576 

2.414 

-.017 

2.576 

2.576 

£.414 

-.017 

-.011 

2.575 

2.414 

2.570 

2.411 

-.017 

-.011 

8x8 

2.413 

-.007 

2.413 

2.412 

-.007 

2.412 

2.406 

-.007 


^Denotes body with singularities in § M (l) the en( ^ 5 ■ 



TABLE III. - DIMENSIONLESS TOTAL ZERO -LIFT WAVE DRAG AT MACH NUMBER 1 FOR SEVERAL CONFIGURATIONS 

OF FIGURE 3 
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TABLE IV.- VALUES OF THE DIMENSIONLESS ZERO-LIFT INTERFERENCE WAVE-DRAG INTEGRAL FUNCTION, l( 
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158942 . 0^5803 .000380 

161698 .050651 .OOO66U 




TABLE IV.- VALUES OF TEE DIMENSIONLESS ZERO-LIFT INTERFERENCE WAVE-DRAG INTEGRAL FUNCTION, l(a- B) - Continued 
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TABLE IV.- VALUES OF THE DIMENSIONLESS ZERO-LIFT INTERFERENCE WAVE-DRAG INTEGRAL FUNCTION, l(a- ■, B) - Continued 
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IV.- VALUES OF THE DIMENSIONLESS ZERO-LIFT INTERFERENCE WAVE-DRAG INTEGRAL FUNCTION, l(a- •, B) - Continued 
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TABLE IV.- VALUES OF THE DIMENSIONLESS ZERO-LIFT INTERFERENCE WAVE-DRAG INTEGRAL FUNCTION, l(a ti , B) - Continued 
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VALUES OF THE DIMENSIONLESS ZERO-LIFT IN^R^RwrE WAVE-DRAG INTEGRAL FUNCTION, iCa^, B) - Concluded 
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Figure 1.- Sketch showing the equal -interval approximation of the second, derivative of the area 

distribution of an equivalent body of revolution. 










Figure 3.- Sketch showing the two basic arrangements of the pairs of bodies of revolution and the 

lengths over which the wave-drag analysis is performed. 
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(b) Half plan forms for 
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